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The Laplace Transform: Introduction 


The Laplace Transform: Excercises 


Continuous-Time Signals: Introduction 
Defines and gives some examples of continuous-time (analog) signals. 


The Merrian-Webster dictionary defines a signal as: 


"A detectable physical quantity or impulse (as a voltage, current, or 
magnetic field strength) by which messages or information can be 
transmitted." 


These are the types of signals which will be of interest in this book. Indeed, 
signals are not only the means by which we perceive the world around us, 
they also enable individuals to communicate with one another on a massive 
scale. So while our primary emphasis in this book will be on the theoretical 
foundations of signal processing, we will also try to give examples of the 
tremendous impact that signals and systems have on society. We will focus 
on two broad classes of signals, discrete-time and continuous-time. We will 
consider discrete-time signals later on in this book. For now, we will focus 
our attention on continuous-time signals. Fortunately, continuous-time 
signals have a very convenient mathematical representation. We represent a 
continuous-time signal as a function x(t) of the real variable t. Here, t 
represents continuous time and we can assign to ¢ any unit of time we deem 
appropriate (seconds, hours, years, etc.). We do not have to make any 
particular assumptions about x(t) such as boundedness (a signal is bounded 
if it has a finite value). Some of the signals we will work with are in fact, 
not bounded (i.e. they take on an infinite value). However most of the 
continuous-time signals we will deal with in the real world are bounded. 
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Temperature signal recorded in Dallas, Texas from Aug. 16 
to Aug. 22, 2002. 


We actually encounter signals every day. Suppose we sketch a graph of the 
temperature outside the Jerry Junkins Electrical Engineering Building on 
the SMU campus as a function of time. The graph might look something 
like in [link]. This is an example of a signal which represents the physical 
quantity temperature as it changes with time during the course of a week. 
[link] shows another common signal, the speech signal. Human speech 
signals are often measured by converting sound (pressure) waves into an 
electrical potential using a microphone. The speech signal therefore 
corresponds to the air pressure measured at the point in space where the 
microphone was located when the speech was recorded. The large 
deviations which the speech signal undergoes corresponds to vowel sounds 
such as “ahhh" or “eeeeh" (voiced sounds) while the smaller portions 
correspond to sounds such as “th" or “sh" (unvoiced sounds). In [link], we 


see yet another signal called an electrocardiogram (EKG). The EKG is a 
voltage which is generated by the heart and measured by subtracting the 
voltage recorded from two points on the human body as seen in [link]. 
Since the heart generates very low-level voltages, the difference signal must 
be amplified by a high-gain amplifier. 
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Measurement of the electrocardiogram 
(EKG). 


Signal Power, Energy, and Frequency 


Signals can be characterized in several different ways. Audio signals 
(music, speech, and really, any kind of sound we can hear) are particularly 
useful because we can use our existing notion of “loudness" and “pitch" 
which we normally associate with an audio signal to develop ways of 
characterizing any kind of signal. In terms of audio signals, we use “power' 
to characterize the loudness of a sound. Audio signals which have greater 
power sound “louder" than signals which have lower power (assuming the 
pitch of the sounds are within the range of human hearing). Of course, 
power is related to the amplitude, or size of the signal. We can develop a 
more precise definition of power. The signal power is defined as: 
Equation: 


' 


1 T/2 
Pz, = lim = / a? (t)dt 


The energy of this signal is similarly defined 


Equation: 
CO 
Cx -| a? (t)dt 


Oo 


We can see that power has units of energy per unit time. Strictly speaking, 
the units for energy depend on the units assigned to the signal. If x(t) is a 
voltage, than the units for e, would be volts?-seconds. Notice also that 
some signals may not have finite energy. As we will see shortly, periodic 
signals do not have finite energy. Signals having a finite energy are 
sometimes called energy signals. Some signals that have infinite energy 
however can have finite power. Such signals are sometimes called power 
signals. 


We use the concept of “frequency” to characterize the pitch of audio 
signals. The frequency of a signal is closely related to the variation of the 
signal with time. Signals which change rapidly with time have higher 


frequencies than signals which are changing slowly with time as seen [link]. 
As we Shall see, signals can also be represented in terms of their 
frequencies, X(j{2), where 2 is a frequency variable. Devices which 
enable us to view the frequency content of a signal in real-time are called 
spectrum analyzers. 
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The signal y(t) contains a greater amount of high frequencies 
than x(t). 


Something to keep in mind is that the signals shown in Figures [link], 
[link], and [link] each have different units (degrees Fahrenheit, pressure, 
and voltage, respectively). So while we can compare relative frequencies 
between these signals, it doesn't make much sense to compare their power 
since each signal has different units. We will take a more formal look at the 
frequency of signals starting in Chapter 2. 


Basic Signal Operations 
Overview of basic continuous-time signal operations. 


We will be considering the following basic operations on signals: 


e Time shifting: 
Equation: 


y(t) = x(t —7) 


The effect that a time shift has on the appearance of a signal is seen in 
[link]. If 7 is a positive number, the time shifted signal, x(t — 7) gets 
shifted to the right, otherwise it gets shifted left. 

e Time reversal: 
Equation: 


Time reversal flips the signal about t = 0 as seen in [link]. 
e Addition: any two signals can be added to form a third signal, 
Equation: 


2(t) = x(t) + y(t) 


e Time scaling: 
Equation: 


y(t) = x(t) 


Time scaling “compresses" the signal if (2 > 1 or “stretches” it if 
2 < 1 (see [link]). 

e Multiplication by a constant, a: 
Equation: 


y(t) = ax(t) 


¢ Multiplication of two signals, their product is also a signal. 
Equation: 


Multiplication of signals has many useful applications in wireless 


communications. 
e Differentiation: 
Equation: 
dx(t 
y(t) = ) 
dt 


e Integration: 
Equation: 


()= [0 


There is another very important signal operation called convolution which 
we will look at in detail in Chapter 3. As we shall see, convolution is a 
combination of several of the above operations. 


(a) x(t) 
7 t 
(b) x(t-T) 


(c) x(-t) 


(a) original signal, (b) time-shift, (c) time-reversal. 


(a) x(t) 


(b) X(Qt) 
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(a) original signal, (b) 2 > 1, (c) 2 < 1. 


Complex Numbers and Complex Arithmetic 


Before we begin studying signals, we need to review some basic aspects of 
complex numbers and complex arithmetic. The rectangular coordinate 
representation of a complex number z is z has the form: 

Equation: 


Z=a+ 7b 


where a and b are real numbers and 7 = /—1. The real part of z is the 
number a, while the imaginary part of z is the number b. We also note that 
7b (jb) = —b? (a real number) since j(7) = —1. Any number having the 
form 

Equation: 


2:= 4b 


where 6 is a real number is an imaginary number. A complex number can 
also be represented in polar coordinates 


Equation: 
z=re” 
where 
Equation: 
t= J a? + b? 


is the magnitude and 
Equation: 


is the phase of the complex number z. The notation for the magnitude and 
phase of a complex number is given by |z| and /z, respectively. Using 
Euler's Identity: 

Equation: 


e*)” —cos (0) + j sin (6) 


it follows that a = r cos (6) and b = r sin (6). [link] illustrates how polar 
coordinates and rectangular coordinates are related. 
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Relationship between 
rectangular and polar 
coordinates. 


Rectangular coordinates and polar coordinates are each useful depending on 
the type of mathematical operation performed on the complex numbers. 
Often, complex numbers are easier to add in rectangular coordinates, but 
multiplication and division is easier in polar coordinates. If z = a+ jbisa 
complex number then its complex conjugate is defined by 

Equation: 


* ° 
z =a-—jb 
in polar coordinates we have 
Equation: 


z =re” 


note that zz = )z|? =r? andz+z = 2a. Also, if 21, z9,..., ZN are 
complex numbers it can be easily shown that 


Equation: 

(21+ 2g t-+++2y) = 2, +2) +-°++2y 
and 
Equation: 


* * * 


(2122°-+ Zn) = 212%9°** ZN 


[link] indicates how two complex numbers combine in terms of addition, 
multiplication, and division when expressed in rectangular and in polar 
coordinates. 


operation rectangular polar 

2 + 29 (a1 + a2) + J (b1 + ba) 

2129 Q1Q_ — b,by + j (aby + aob,) rpreed (+9) 

21/29 Sie Ah 1 @i(O—%) 
2+» 


Operations on two complex numbers, z; = a, + 7b; = r,e7 and 
Zy = ay + jbo = rye)”, The sum of two complex numbers is cumbersome 
to express in polar coordinates, and is not shown. 


Periodic Signals 


Periodic signals have the following property: 
Equation: 


a(t) = «(t+ kT) 


where & is an integer and 7’ is called the fundamental period. Periodic have 
the property that they “repeat" every 7’ seconds. For periodic signals, the 
power can be defined as 

Equation: 


[link] shows an example of a periodic signal. We will study the frequency 
content of periodic signals in some detail in Chapter 2. 


General periodic signal. 


Sinusoidal Signals 


Sinusoidal Signals 


Sinusoidal signals are perhaps the most important type of signal that we 
will encounter in signal processing. There are two basic types of signals, the 
cosine: 

Equation: 


a(t) = A cos (Mt) 


and the sine: 
Equation: 


z(t) = A sin (Mt) 


where A is a real constant. Plots of the sine and cosine signals are shown in 
[link]. Sinusoidal signals are periodic signals. The period of the cosine and 
sine signals shown above is given by T’ = 27/2. The frequency of the 
signals is (2 = 27/T which has units of rad/sec . Equivalently, the 
frequency can be expressed as 1/7’, which has units of sec, cycles/sec, or 
Hz. The quantity (2¢ has units of radians and is often called the phase of the 
sinusoid. Recalling the effect of a time shift on the appearance of a signal, 
we can observe from [link] that the sine signal is obtained by shifting the 
cosine signal by T’/4 seconds, i.e. 
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Cosine and sine signals. Each signal is periodic with 
period T = 27/9. 


Equation: 

sin (2t) =cos (Q(t — T/4)) 
and since T’ = 27/2, we have 
Equation: 

sin (2t) =cos (Mt — 7/2)) 
Similarly, we have 
Equation: 


cos (2t) =sin (2 + 7/2)) 


Using Euler's Identity, we can also write: 


Equation: 
A (5008 — jt 
A cos (t) = ml + eM) 
and 
Equation: 
: A (5001 — jt 
A sin (t) = — (e™ — e ) 
2j 
The quantity e7 is called a complex sinusoid and can be expressed as 
Equation: 


et IM —cos (Qt) + 7 sin (Mt) 


There are a number of trigonometric identities which are sometimes useful. 
These are shown in [link]. [link] shows some basic calculus operations on 
sine and cosine signals. 


sin (0) =cos (0 — 77/2) 
cos (9) =sin (0 + 7/2) 
sin (0,) sin (0,) = +[cos (6; — 42)— cos (0; + 4)| 
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sin (8,) cos (62) = $ [sin (0; — 6,)— sin (0; + 4)| 


cos (0;) cos (02) = 1 [cos (0; — 02)+ cos (8; + 4.)| 
a cos (9) + bsin (8) = Va? + b?cos(6— tan! (4)) 
cos (0; + 02) =cos (01) cos (82) sin (01) sin (62) 
sin (0; + 62) =sin (01) cos (02)+ sin (61) cos (02) 


Useful trigonometric identities. 


— cos (2M) = —N sin (Mt) 
— sin (2t) = 2 cos (Nt) 
f cos (Qt)dt = = sin (2t) 
Jf sin (Qt)dt = —~, cos (Mt) 
J, sin (kQot) cos (nQot)dt = 0 
Jo. sin (kQot) sin (nAot)dt =0,k An 
Jo. cos (kQot) cos (nQot)dt = 0,k An 
i sin? (nQ,t)dt = T/2 
i cos? (nQ,t)dt = T/2 


Derivatives and integrals of sinusoidal signals. 


Now suppose that we have a sum of two sinusoids, say 
Equation: 


x (t) =cos (2,t)+ cos (ot) 


It is of interest to know what the period 7’ of the sum of 2 sinusoids is. We 
must have 
Equation: 


a(t—T) = cos (M(t —T))+ cos (MQ (t — T)) 
cos (2)t — 2;T)+ cos (Qt — 22T) 


It follows that 2,;T = 27k and 22T = 27, where k and | are integers. 
Solving these two equations for T’ gives T = 27k/Q) = 2n1/Q2. We wish 
to select the shortest possible period, since any integer multiple of the 
period is also a period. To do this we note that since 27k/Q, = 27l/o, 
we can write 

Equation: 


2} k 


Q 1 


so we seek the smallest integers & and / that satisfy [link]. This can be done 
by finding the greatest common divisor between k and /. For example if 

92, = 107 and (22g = 157, we have k = 2 and / = 3, after dividing out 5, 
the greatest common divisor between 10 and 15. So the period is 

T = 2nk/Q, = 0.4 sec. On the other hand, if 2; = 107 and 2) = 10.17 
, we find that k = 100 and / = 101 and the period increases to 

T = 2nk/Q, = 20 sec. Notice also that if the ratio of 2, and 2 is not a 
rational number, then x(t) is not periodic! 


If there are more than two sinusoids, it is probably easiest to find the period 
of one pair of sinusoids at a time, using the two lowest frequencies (which 
will have a longer period). Once the frequency of the first two sinusoids has 


been found, replace them with a single sinusoid at the composite frequency 
corresponding to the first two sinusoids and compare it with the third 
sinusoid, and so on. 


Introduction to Continuous-Time Signals: Exercises 


1. Consider the signals shown in [link]. Sketch the following 


signals|[ footnote]: 
a. £1 (t) + 22 (Et). 
b. Ly (— t) — @ (t — 1). 
c. x, (-t+1). 
d. “2 (t = 1); 
e. £1 (2t). 
f. Ly (t/2). 
g. x1 (t)x3 (t+ 1). 
h. L3 (2t = A). 
i. x3 (— 2t — A). 


Assume that for step discontinuities, the signal takes on the greater of 
the two values. 
2. For each of the signals in [link]: 


a. What is the period? 
b. Sketch x(t — 0.25), and x(t + 1). 
c. Find the power for each signal. 


3. Suppose that z3 = 3+ 92, 2 = 4+ 75. Find: 


Z1|, ys |Z2l, L290 

z, + 22 in rectangular coordinates. 

2122 in rectangular and polar coordinates. 
21/2Z2 in polar coordinates. 


o Oo 0 °O 


4. What is the period, frequency, and power of the sinusoidal signal 
elt) = 2 cos (50). 

5. Can you find a general formula for the power of the sinusoidal signal 
x(t) = A cos (M2)? 

6. Express 2 cos (10t) + 3 sin (10¢) as a single sinusoidal signal. 

Sketch x(t) =sin (¢ — 6) for 0 = 1/4, 7/2, and 37/4. 

8. Sketch x(t) =sin (2t¢ — 0) for 6 = 1/4, 2/2, and 37/4. 


= 


9. Consider the motion of the second hand of a clock. Assume the length 
of the second hand is 1 meter. (a) what is the angular frequency of the 
second hand. (b) find an expression for the horizontal and vertical 
displacements of the tip of the second hand, assuming the origin is at 
the clock center and t = 0 when the second hand is over the 3. 
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Signals for problem 2. 


Equation: 


Symmetry Properties of Periodic Signals 


A signal has even symmetry of it satisfies: 
Equation: 


and odd symmetry if it satisfies 
Equation: 


x(t) = —ax(-t) 


[link] shows pictures of periodic even and odd symmetric signals. If x(t) is 
an odd symmetric periodic signal, then we must have: 


Equation: 
to+T 
i x (t)dt = 0 
t 


0 


This is easy to see if we choose tp = —T/2. 


(b) 


(a) Even-symmetric, and (b) odd-symmetric 
periodic signals. Note that the integral over any 
period of an odd-symmetric periodic signal is zero. 


We also note that the product of two even signals is also even while the 
product of an even signal and an odd signal must be odd. Finally, the 
product of two odd signals must be even. For example, suppose 2, (t) has 
odd symmetry and x, (t) has even symmetry. Their product has odd 
symmetry because if y(t) = x, (t)x, (t), then 

y(—t) = &o(—t)&e (—t) = —y(t). 


Trigonometric Form of the Fourier Series 


A major goal of this book is to develop tools which will enable us to study 
the frequency content of signals. An important first step is the Fourier 
Series. The Fourier Series enables us to completely characterize the 
frequency content of a periodic signal[footnote]. A periodic signal x(t) can 
be expressed in terms of the Fourier Series, which is given by: 

There are periodic signals for which a Fourier series doesn't exist, 
conditions for existence of the Fourier series are given below. 

Equation: 


z(t) =ao+ S> An cos (Not) + S- b, sin (nMot) 


where 
Equation: 


{0 = 


is the fundamental frequency of the periodic signal. Examination of [link] 
suggests that periodic signals can be represented as a sum of suitably scaled 
cosine and sine waveforms at frequencies of 29, 2929, 3p, .... The cosine 
and sine terms at frequency np are called the n™ harmonics. Evidently, 
periodic signals contain only the fundamental frequency and its harmonics. 
A periodic signal cannot contain a frequency that is not an integer multiple 
of its fundamental frequency. 


In order to find the Fourier Series, we must compute the Fourier Series 
coefficients. These are given by 
Equation: 


Equation: 


) to+T 
a, = = x(t) cos (nQot)dt, n=1,2,... 
T Jt, 
Equation: 
) to+T 
bp, = — x(t) sin (nQot)dt, n=1,2,... 
Tle 


From our discussion of even and odd symmetric signals, it is clear that if 
z(t) is even, then x (t) sin (n{2ot) must be odd and so b,, = 0. Also if, 
z(t) has odd symmetry, then x (¢) cos (not) also has odd symmetry and 
hence a, = 0 (see exercise [link]). Moreover, if a signal is even, since 

x (t) cos (n{ot) is also even, if we use the fact that for any even 
symmetric periodic signal u(t), 


Equation: 

T/2 T/2 

i v(t)dt = 2 / v (t)dt 
_T/2 0 


then setting tg = —T'/2 in [link] gives, 
Equation: 


4 T/2 
tn = =f x(t) cos (nQot)dt, n=1,2,... 
0 


This can sometimes lead to a savings in the number of integrals that must 
be computed. Similarly, if 2(¢) has odd symmetry, we have 
Equation: 


4 T/2 
b, = =| z(t) sin (nQot)dt, n=1,2,... 
0 


Example 2.1 Consider the signal in [link]. This signal has even symmetry, 
hence all of the 6, = 0. We compute apo using, 
Equation: 


1 totT 
ag = — x (t)dt 
LT Jt, 


which we recognize as the area of one period, divided by the period. Hence, 
ao = T/T. Next, using [link] we get 
Equation: 


» q]2 


Qn = 7 - cos (nQot)dt 


Note how the limits of integration only go from —7/2 to 7/2 since x(t) is 
zero everywhere else. Evaluating this integral leads to 
Equation: 


_ 27 sin(nQrT/2) 


= =1,2,... 
eT alQyr/2 7° 


[link] shows the first few Fourier Series coefficients for 7 = 1/2 and 

T = 1. If we attempt to reconstruct x(t) based on only a limited number 
(say, V) of Fourier Series coefficients, we have 

Equation: 


x(t) =ao+ 


N 
An Cos (Not) 


=0 


Figures [link] and [link] show & (t) for VN = 10, and N = 50, respectively. 
The ringing characteristic is known as Gibb's phenomenon and disappears 
only as N approaches oo. 


The following example looks at the Fourier series of an odd-symmetric 
signal, a sawtooth signal. 


Example 2.2 Now let's compute the Fourier series for the signal in [link]. 
The signal is odd-symmetric, so all of the a, are zero. The period is 

T = 3/2, hence 2 = 47/3. Using [link], the 6, coefficients are found by 
computing the following integral, 

Equation: 


g pi 
a sf t sin (4rnt /3)dt 
3 J_1/2 


After integrating by parts, we get 
Equation: 


i 2 
— 3 sin uae), _ 9 £08 (27n/3) | dreaded. 
(7n) TM 


These are plotted in [link] and approximations of x(t) using N = 10 and 
N = 50 coefficients are shown in Figures [link] and [link], respectively. 
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Half-Wave Symmetry 


Half-Wave Symmetry 


Periodic signals having half-wave symmetry have the property 
Equation: 


x(t) —2(t — T/2) 
a(t) = —a(t+T7/2) 


It turns out that signals with this type of symmetry only have odd-numbered 
harmonics, the even harmonics are zero. To see this, lets look at the formula 
for the coefficients a,,: 


Equation: 
) tot+T 
an = x (t) cos (nQot)dt 
to 
) to+T/2 tot+T 
- — x (t) cos (nQot)dt + / x (t) cos (nQot)dt 
T to tot+T/2 
2 
= ql = Ip| 


Making the substitution 7 = ¢ — T’/2 in Iz gives 
Equation: 


totT'/2 
e-= / GID coe Gs PRD de 


0 


to+T/2 
a -| x (T) cos (nM (7 + T/2))dr 


0 


The quantity cos (no (7 + T'/2)) =cos (nr + nz) can be simplified 
using the trigonometric identity 
Equation: 


cos (u + v) =cos (u) cos (v)+ sin (u) sin (v) 
We have 
Equation: 


cos (nQ7 + nm) =cos (n2T) cos (n7)— sin (nT) sin (nz) 
= (—1)” cos (nr) — 0 


Therefore 
Equation: 


Ip x (rT) cos (nQ2Qr)dr 


| 
| 
| 
——— 
3 
=> 
S 
+ 
y 
SS 
i) 


and we can write: 
Equation: 


to tT /2 
i= = —(-1)") i mice ital 


0 


From this expression we find that a,, = 0 whenever n is even. In fact, we 
have 
Equation: 


to 
0, n, even 


One i uae z (t) cos (nQot)dt, n, odd 


A similar derivation leads to 


Equation: 


T Jto 
0, n, even 


b -{3 ht” g(t) sin (nQot)dt, n, odd 


A good choice of to can lead to a considerable savings in time when 
calculating the Fourier Series of half-wave symmetric signals. Note that 
half-wave symmetric signals need not have odd or even symmetry for the 
above formulae to apply. If a signal has half-wave symmetry and in addition 
has odd or even symmetry, then some additional simplification is possible. 
Consider the case when a half-wave symmetric signal also has even 
symmetry. Then clearly b,, = 0, and [link] applies. However since the 
integrand in [link] is the product of two even signals, x(t) and cos (n{ot), 
it too has even symmetry. Therefore, instead of integrating from, say, —T’/4 
to T’/4, we need only integrate from 0 to T’/4 and multiply the result by 2. 
Therefore the formula for a, for an even, half-wave symmetric signal 
becomes: 

Equation: 


a, = i Pome: (t) cos (nQt)dt, n, odd 


0, n, even 
Equation: 
b, = 0 
For an odd half-wave symmetric signals, a similar argument leads to 
Equation: 
i, =—U 


Equation: 
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J ie m 


(t) sin (n{Qot)dt, n, odd 


nm, even 


Convergence of the Fourier Series 


Consider the trigonometric form of the Fourier series 
Equation: 


z(t) =ag+ > An cos (Not) + S- b, sin (nMot) 


n=1 n=1 


It is important to state under what conditions this series (the right-hand side 
of [link]) will actually converge to x(t). The nature of the convergence also 
needs to be specified. There are several ways of defining the convergence 
of a series. 


1. Uniform convergence: define the finite sum: 


Equation: 
N N 
ty (t)=ao+ > ay cos (nMot) + >: b,, sin (nt) 
n=1 n=1 


where JN is finite. Then the series converges uniformly if the absolute 
value of x (t) — xy (t) satisfies 
Equation: 


|x (t) — an (t)| <e€ 


for all values of t and some small positive constant e. 
2. Point-wise convergence: as with uniform convergence, we require that 
Equation: 


|x (t) — LN(t) (t)| <— 


for all ¢. The main difference between uniform and point-wise 
convergence is that for the latter, the number of terms in the 
summation V(t) needed to get the error below € may vary for different 
values of t. 

3. Mean-squared convergence: here, the series converges if for all ¢: 


Equation: 


tot+T 9 
lim |x (t) — nq (¢)| dt = 0 


N-oo t 
0 


Gibb's phenomenon, mentioned in some of the examples above, is an 
example of mean-squared convergence of the series. The overshoot in 
Gibb's phenomenon occurs only at abrupt discontinuities. Moreover, 
the height of the overshoot stays the same independently of the number 
of terms in the series, VV. The overshoot merely becomes less 
noticeable because it becomes more and more narrow as JN increases. 


Dirichlet has given a series of conditions which are necessary for a periodic 
signal to have a Fourier series. If these conditions are met, then 


¢ the Fourier series has point-wise convergence for all t at which x(t) is 
continuous. 

e where x(t) has a discontinuity, then the series converges to the 
midpoint between the two values on either side of the discontinuity. 


The Dirichlet Conditions are: 


1. x(t) has to be absolutely integrable on any period: 


Equation: 
to+T 
/ a(i\|dé < 60 
to 


2. a(t) can have only a finite number of discontinuities on any period. 
3. 2(t) can have only a finite number of extrema on any period. 


Most periodic signals of practical interest satisfy these conditions. 
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Complex Form of the Fourier Series 


The trigonometric form of the Fourier Series, shown in [link] can be 
converted into a more convenient form by doing the following 
substitutions: 

Equation: 

einQot 4 e InMot 


cos (not) = ; 
Equation: 
eInMot = e InMot 
in (not) = —————_ 
sin (n{Mot) 2 


After some straight-forward rearranging, we obtain 
Equation: 


n—J qb n bn 
z(t) =a0+ | oP “a emo 4 yo) One te |e —jnQot 


Keeping in mind that a,, and by, are only defined for positive values of n, 
lets sum over the negative integers in the second summation: 
Equation: 


o(t)=a0+ [5 ms emmy S| a 2 | em 


n=-1 


Next, let's assume that a,, and 6, are defined for both positive and negative 
n. In this case, we find that a,, = a_, and b,, = —b_,, since 
Equation: 


) tot+T 
an = a. x (t) cos (—not)dt 


totT 
= — / x (t) cos (nQot)dt 
t 


and 
Equation: 


tot+T 
bb», = zi x (t) sin (—n{2ot)dt 
T Jet, 


) totT 
= <i x (t) sin (nQt)dt 
to 
=e 


Using this fact, we can rewrite [link] as 
Equation: 


x(t) = a+ [25% | eine S ae *|e gn Qot 


n=—-1 


If we define[footnote|cg = ag, and 


The notation “=" is often used instead of the “=" sign when defining a new 
variable. 
Equation: 
c, = an — Jon 
_— 2 


then we can rewrite [link] as 
Equation: 


CO 
x(t) = S- c, esreot 


n=—OCoO 


which is called the complex form of the Fourier Series. Note that since 
QA_n — Q, and b_, — —b,, we have 
Equation: 


This means that 


Equation: 
\C_n| = |en| 
and 
Equation: 
LC = —ZCp 


Next, we must find formulas for finding the c, given x(t). We first look at 
a property of complex exponentials: 


Equation: 
eel = 
| ot eth Qot ge — Tk) 
to 0, otherwise 


To see this, we note that 
Equation: 


to+T tot+T tot+T 
/ ero ay = / cos (kQot) + ii, sin (kQot)dt 


to to to 


It's easy to see that k{29 also has period 7’, hence the integral is over k 
periods of cos (k{2ot) and sin (kot). Therefore, if k ~ 0, then 


Equation: 
tot+T 
/ Shot dt — 0) 
to 


otherwise 
Equation: 


totT ; to tT 
/ eIhOot qe — / dt —T 
to to 


We use [link] to derive an equation for c, as follows. Consider the integral 


Equation: 
totT ; 
/ a (t)e 9 dt 
t 


0 


Substituting the complex form of the Fourier Series of x(¢) in [link], (using 
k as the index of summation) we obtain 


Equation: 
tot+T (oe) 
i > Cpe dhrrt e InQot ay 
t 


0 k=—00 


Rearranging the order of integration and summation, combining the 
exponents, and using [link] gives 
Equation: 


[o-e) totT 
S- Cn / ellk-n)Qot dt — Te, 
to 


k=—0o 


Using this result, we find that 
Equation: 


1 totT 
Ci == a (t)e 9" de 


Example 2.1 Let's now find the complex form of the Fourier Series for the 
signal in Example [link]. The integral to be evaluated is 
Equation: 


9 705 = 
——— = / Qte I3™ dt 
3 J_o.5 


Integrating by parts yields 
Equation: 


= + cos (27n/3) — 5 sin (27/3) 
nt 


(nm) 


[link] shows the magnitude of the coefficients, |c,,|. Note that the complex 
Fourier Series coefficients have even symmetry as was mentioned earlier. 


Fourier Series coefficients for Example "Complex Form of 


the Fourier Series". 


Can the basic formula for computing the c, in [link] be simplified when 
z(t) has either even, odd, or half-wave symmetry? The answer is yes. We 
simply use the fact that 

Equation: 


An — jon 
Ch = —— 
2 


and solve for a, and 6,, using the formulae given above for even, odd, or 
half-wave symmetric signals. This avoids having to integrate complex 
quantities. This can also be seen by noting that (setting tg = T'/2 in [link]): 
Equation: 


CG =] = x (t)e Irot dt 

ie eee (t) 
1 pri T/2 

= = x (t) cos (nQMt)dt — Zi ) sin (nQot)dt 
PT J_ryo ri. 
1 j 

= Zan — —bn 

5 


Alternately, if 7(¢) has half-wave symmetry, we can use [link], [link], and 
[link] to get 
Equation: 


a= {? 2 frig (t) e Mmtdt mn odd 
0, 


nm even 


Unlike the trigonometric form, we cannot simplify this further if x(t) is 
even or odd symmetric since e~ 4” has neither even nor odd symmetry. 


Example 2.2 In this example we will look at the effect of adjusting the 
period of a pulse train signal. Consider the signal depicted in [Link]. 


x(t) 


a i -/2 0 7/2 fi 


Pulse train having period 7’ used in Example "Complex Form of the 
Fourier Series". 


The Fourier Series coefficients for this signal are given by 


Equation: 


Cn 


[link] shows the magnitude of |c, 


1 F/2 
acs e InLot dt 
si —7/2 
4 . 
—jnQt/2 ane) 
jnQoT (¢ 
7 sin (nQ%r7/2) 
T Not /2 
arsine (nQor /2) 


, the amplitude spectrum, for T’ = 1 and 


T = 1/2. as well as the Fourier Series for the signal based on the first 30 


coefficients 
Equation: 


30 
Z(t) = S- jem 


n=—30 


Similar plots are shown in Figures [link], and [link], for 7’ = 4, and 7’ = 8, 


respectively. 


Example "Complex Form of the Fourier Series", T = 1,7 = 1/2: 
(top) Fourier Series coefficient magnitudes, (b) Z (t). 


()x 


t (sec) 


Example "Complex Form of the Fourier Series", T = 4,7 = 1/2: 
(top) Fourier Series coefficient magnitudes, (b) Z (t). 


x(t) 


Example "Complex Form of the Fourier Series", T = 8,7 = 1/2: 
(top) Fourier Series coefficient magnitudes, (b) % (t). 


0.07 


x(t) 


0.4 


0.2 


This example illustrates several important points about the Fourier Series: 
As the period T increases, {29 decreases in magnitude (this is obvious since 
Qo = 27/T). Therefore, as the period increases, successive Fourier Series 
coefficients represent more closely spaced frequencies. The frequencies 
corresponding to each n are given by the following table: 


n 22 

0 0 

=a | +2 
2 +2929 
+n +n 


This table establishes a relation between n and the frequency variable §2. In 
particular, if 7’ = 1, we have (29 = 27 and 


+1 +27 


az +47 


—n +2n7 


If T = T, then 2Q9 = 7/2 and 


n 2 

0 0 

+1 +n/2 
ea aay |) 
+n +n /2 


and if T = 8, we have 29 = 7/4 and 


+1 +n/4 


+2 +n/2 


+n +n7/4 


Note that in all three cases, the first zero coefficient corresponds to the 
value of n for which (2 = 47. Also, as T gets bigger, the c,, appear to 
resemble more closely spaced samples of a continuous function of 
frequency (since the n{2 are more closely spaced). Can you determine what 
this function is? 


As we Shall see, by letting the period T’ get large (infinitely large), we will 
derive the Fourier Transform in the next chapter. 
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Parseval's Theorem for the Fourier Series 


Recall that in Chapter 1, we defined the power of a periodic signal as 
Equation: 


1 totT 


2 
iu: = > x (t)dt 
LT Ji, 


where 7’ is the period. Using the complex form of the Fourier series, we can 
write 
Equation: 


a(t)? = ( S- Tad ( 3 ad 


n=—CcoO M=— CoO 


where we have used the fact that z(t)” = x (t)x(t) , Le. since x(t) is real 
a (t) = x(t) . Applying [link] and [link] gives 


Equation: 
OO 00 ‘ 
( S° aco) ( S° cere) 
n=—0o m=—oo 
OO  . 
S° S° Cy Cp, Cm) Mot 


n=—-W M=—CcoO 


oo : ae 
S } Cr +) "CnC er mot 


n=— oo nAém 


a(t)? 


Substituting this quantity into [link] gives 
Equation: 


to 


Pz: = 
n=—0oo nz=m 
(oe) 1 to+T re ; 
- Stated [Secon 
n=—0o - to ném 


It is straight-forward to show that 
Equation: 


totT . 
Ew i S° CC Coe "dt = 0 
L to nz=m 


This leads to Parseval's Theorem for the Fourier series: 
Equation: 


00 
Pz = > le,| 


n=— Oo 


which states that the power of a periodic signal is the sum of the magnitude 
of the complex Fourier series coefficients. 


The Fourier Series: Exercises 


1 


Z; 


3 


. Show that an even-symmetric periodic signal has Fourier Series 


coefficients b,, = 0 while an odd-symmetric signal has a, = 0. 

Find the trigonometric form of the Fourier Series of the periodic signal 
shown in [link]. 

Find the trigonometric form of the Fourier Series for the periodic 
signal shown in [link]. 


. Find the trigonometric form of the Fourier Series for the periodic 


signal shown in [link] for 7 = 1, JT’ = 10. 


. Suppose that z(t) = 5 + 3 cos (5t) — 2 sin (3t)+ cos (45t). 


1. Find the period of this periodic signal. 
2. Find the trigonometric form of the Fourier Series. 


. Find the complex form of the Fourier Series of the periodic signal 


shown in [link]. 


. Find the complex form of the Fourier Series of the periodic signal 


shown in [link]. 


. Find the complex form of the Fourier Series for the signal in [link] 


using: 


Lares 1 
ie ea ad 
For each case plot the magnitude of the Fourier Series coefficients. 
You may use Matlab or some other programming language to do this. 


. Show that 


Equation: 


1 tot+T * 
— / S° CnC Coro dt = 0 
‘ 


0 ném 


x(t) 


x(t) 


A 0 4 5 T 


Pulse train signal for problems [link] and [link]. 


Derivation of the Fourier Transform 
The Fourier transform is derived from the Fourier series. 


Let's begin by writing down the formula for the complex form of the 
Fourier Series: 
Equation: 


CO 
x(t) = i cperreot 


n=— Oo 


as well as the corresponding Fourier Series coefficients: 
Equation: 


1 totT ; 
=> a (t)e 30 dt 
Tha 


As was mentioned in Chapter 2, as the period 7’ gets large, the Fourier 
Series coefficients represent more closely spaced frequencies. Lets take the 
limit as the period T goes to infinity. We first note that the fundamental 
frequency approaches a differential 


Equation: 
27 
26 = df? 
Fe 
consequently 
Equation: 
1 4 di) 


— = 
di 2 21 


The nth harmonic, n{2o, in the limit approaches the frequency variable 2 
Equation: 


nQy > 2 


From equation [link], we have 
Equation: 


cP f a (t)e dt 


The right hand side of [link] is called the Fourier Transform of x(t): 
Equation: 


X (jQ2) = i a (t)e 1 dt 


[o,@) 


Now, using [link], [link], and [link] in equation [link] gives 
Equation: 


1 
a / X (jean 
OR J 56 


which corresponds to the inverse Fourier Transform. Equations [link] and 
[link] represent what is known as a transform pair. The following notation 
is used to denote a Fourier Transform pair 

Equation: 


a(t) o X(j2) 


We say that x(t) is a time domain signal while X(j22) is a frequency 
domain signal. Some additional notation which is sometimes used is 
Equation: 


X (j2) = Fi{a(t)} 


and 
Equation: 


x (t) = F{X(jQ)} 
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Properties of the Fourier Transform 
Several of the most important properties of the Fourier transform are derived. 


Properties of the Fourier Transform 
The Fourier Transform (FT) has several important properties which will be useful: 


1. Linearity: 
Equation: 


ax, (t) + Bre (t) 6 aX, (J) + BX2 (j22) 


where q@ and are constants. This property is easy to verify by plugging the 
left side of [link] into the definition of the FT. 

. Time shift: 
Equation: 


N 


a(t—T) ve 9°" X (52) 


To derive this property we simply take the FT of x(t — rT) 
Equation: 


/ a(t—r)e Idt 


CO 


using the variable substitution y = t — 7 leads to 


Equation: 
CY opr 
and 
Equation: 
ay=di 


We also note that if £ = +oo then 7 = too. Substituting [link], [link], and 
the limits of integration into [link] gives 
Equation: 


- a(yje wma, = / x (ye IM dy 


[e@) CO 


e I°7X (52) 


which is the desired result. 
. Frequency shift: 
Equation: 


x (t)e¥" <> Xj (2 — N)) 


Deriving the frequency shift property is a bit easier than the time shift 
property. Again, using the definition of FT we get: 
Equation: 


/ a (t)eI eI dt - | a (t)e 1PM ae 


[e@) CO 


X(j(2— M)) 


. Time reversal: 
Equation: 


a(—t) 4+ X(—j2) 


To derive this property, we again begin with the definition of FT: 


Equation: 
Co # 
/ a (—t)e dt 


[e.@) 


and make the substitution y = —t. We observe that dt = —dy and that if the 
limits of integration for t are -Eoo, then the limits of integration for yy are Fy 
. Making these substitutions into [link] gives 

Equation: 


7 / x(yJerdy = / x (yer "dy 


= X(-J0) 


Note that if x(t) is real, then X (—jQ) = X(jQ) . 


5. Time scaling: Suppose we have y(t) = z(at),a > 0. We have 
Equation: 


Go) = a a (at)e "dt 


CO 


Using the substitution y = at leads to 
Equation: 


Y(j2) 


= of (yJe~ 52/0 gy 


(2 


6. Convolution: The convolution integral is given by 
Equation: 


git) = io x (r)h(t —7r)dr 


CO 


The convolution property is given by 
Equation: 


Y(jQ) > XGQ)H(G2) 


To derive this important property, we again use the FT definition: 
Equation: 
oO . 
[. y (the dt 


[. [= (r)h(t —r)e 9 drdt 
[2@|f_ne-newmae dr 


Using the time shift property, the quantity in the brackets is e 907 (722), 
giving 
Equation: 


Y(j2) 


¥(jQ) = [. * (re JO" (jQ)dr 


= (50) [ 2( (r)e I?" dr 
H(gj2)X (92) 


Therefore, convolution in the time domain corresponds to multiplication in 
the frequency domain. 

7. Multiplication (Modulation): 
Equation: 


w(t) =2@y) > = [ xG(2-e)¥ Gea 


Notice that multiplication in the time domain corresponds to convolution in 
the frequency domain. This property can be understood by applying the 
inverse Fourier Transform [link] to the right side of [link] 

Equation: 


w(t) =f .3 “2 fx X (j(Q— @))Y (jO@)e?dodN 


- =f. Y ( 0) | 5 [| xG(a-eyeao dO 


The quantity inside the brackets is the inverse Fourier Transform of a 
frequency shifted Fourier Transform, 
Equation: 


1 me 
wii) = a / Y (jO) |x (t)e?"] dO 
7 —=00 
1 a 
= a(t)— | Y(jO)e%"do 
cs f _YUe)e 
= a(t)y(t) 
8. Duality: The duality property allows us to find the Fourier transform of time- 
domain signals whose functional forms correspond to known Fourier 


transforms, X (jt). To derive the property, we start with the inverse Fourier 
transform: 


Equation: 
a(t) = aa / “ X (jQ)ei*dA 
2M J es 


Changing the sign of ¢ and rearranging, 
Equation: 


Qra (—t) = / X (jQ)e I" dQ 


Now if we swap the ¢ and the £2 in [link], we arrive at the desired result 
Equation: 


Qne (—2) = / X (jt)e iat 


The right-hand side of [link] is recognized as the FT of X(jt), so we have 
Equation: 


X(jt) @ 2n2(—2) 


The properties associated with the Fourier Transform are summarized in [link]. 


Property y(t) Y(jQ) 

Linearity ax, (t) + Bre (t) aX, (72) + BX2 (j22) 
Time Shift a(t — T) X (jQ)e I" 

F : 

Frequency a (t)ei*! X(j(@—M)) 

Time z(—t) X (—j2) 


Reversal 


ae (at) qe) 


Scaling a 

Convolution x(t)*h(t) X (j2)A (72) 

Modulation z(t)w(t) = [= X (7 (2 — @))W (j@)dO 
Duality X (jt) 2rx(—22) 


Fourier Transform properties. 


Symmetry Properties of the Fourier Transform 
Looks at some of the consequences of conjugate symmetry of the Fourier 
transform of real signals. 


When 2(t) is real, the Fourier transform has conjugate symmetry, 
X (—jQ) = X(jQ) . It is not hard to see this: 
Equation: 


* 


X(jQ) = / _ 2 (eat 


Oo 


= ‘ * [x (te 3] “at 


OO 


[- a (t)er dt 
X(—ji) 


where the second equality uses the definition of a Riemann integral as the 
limiting case of a summation, and the fact that the complex conjugate of a 
sum is equal to the sum of the complex conjugates. The third equality used 
the fact that the complex conjugate of a product is equal to the product of 
complex conjugates. 


Letting X (72) = a(j) + jb(7M2), it follows that 
Equation: 


X(—j2) = a(—jQ) + jb(—Jj2) 


and 
Equation: 


X(j2) = a (92) — jb (52) 


Equating [link] and [link] gives a(jQ) = a(—j2) and b(—j2) = —b(j2), 
which implies that the real and imaginary parts of X(j2) have even and 
odd symmetry, respectively. A consequence of this is that 

AG) = X(j2)"| = |X(—Jj.)|, that is, the magnitude of the Fourier 
transform has even symmetry. It can similarly be shown that the phase of 
the Fourier transform has odd symmetry. 


The Unit Impulse Function 


The Unit Impulse Function 


The unit impulse is very useful in the analysis of signals, linear systems, 
and sampling. Consider the plot of a rectangular pulse in [link]. Note the 
height of the pulse is 1/7 and the width of the pulse is 7. So we can write 


Equation: 
[e,@) 
/ zp (t)dt = 1 


CO 


As we let T get small, then the width of the pulse gets successively 
narrower and its height gets progressively higher. In the limit as 7 
approaches zero, we have a pulse which has infinite height, and zero width, 
yet its area is still one. We define the unit impulse function as 
Equation: 

6 (t) =lim z, (t) 


tT 0 


x,(t) o(t) 
L/t 
T— 0 
> 
t t 
-/2 0 12 0 


Rectangular pulse, x, (t) approaches the unit impulse 
function, 6(t), as 7 approaches zero. 


The area under 6(t) is one, and so we can write 
Equation: 


[ se-nae=1 


If we multiply the unit impulse by a constant, K, its area is now equal to 
that constant, i.e. 
Equation: 


/ Ké(t—1)dt =K 


The area of the unit impulse is usually indicated by the number shown next 
to the arrow as seen in [link]. 


T 


Ké(t — 7). 


Suppose we multiply the signal x(t) with a time-shifted unit impulse, 
6(¢ — r). The product is a unit impulse, having an area of x(r). This is 
illustrated in [link]. 


O(t - T) 


t 
T 
X 
x(t) 
t 
lI 
x(T) 
t 
T 


Sifting property of unit 
impulse, the product of 
the two signals, x(t) 
and 6(t — 7), is 
x(r)d(t — 7). 


Consequently, the area 


under x(7)d(t — T) is 


ar): 


In other words, 
Equation: 


Equation [link] is called the sifting property of the unit impulse. As we will 
see, the sifting property of the unit impulse will be very useful. 


The Unit Step Function 


The unit step function is defined as 
Equation: 


This function is useful for defining signals which we wish to start at t 
In other words, often, we would like for signals to be zero for negative 
values of t. We can force this situation by simply multiplying by u ¢ . 


Fourier Transform of Common Signals 
Fourier transform of common signals are derived. 


Fourier Transform of Common Signals 


Next, we'll derive the FT of some basic continuous-time signals. [link] 
Summarizes these transform pairs. 


Rectangular pulse 


Let's begin with the rectangular pulse 
Equation: 


lL, $27 /2 
0, t>7/2 


rect (t,7) = { 


The pulse function, rect(t, 7) is shown in [link]. Evaluating the Fourier 
transform integral with x(t) = rect(t, 7) gives 
Equation: 


7/2 
X(jQ2) = 7 eI" dt 
—7/2 
— a1 oi ie 
982 
_ 2 ew _ eal 
982 
sin (Q7/2) 
Qr /2 
= rsinc(27/2) 


—7/2 


A plot of X(j) is shown in [link]. 


x(t) = rect(t, T) 


1 


ies i “or an [2m 0 2m\ /an OR 
T OT T t t 26 


Fourier transform pair showing the rectangular pulse signal (left) and 
its Fourier transform, the sinc function (right). 


Note that when 2 = 0, X(j2) = 7. We now have the following transform 
pair: 

Equation: 

sin (27/2) 


rect (t,T) OT Or)? 


Impulse 


The unit impulse function was described in a previous section. From the 
sifting property of the impulse function we find that 
Equation: 


X(jQ2) = / ” § (te iat 
= 1 


or 
Equation: 


Complex Exponential 


The complex exponential function, x (t) = e7¢, has a Fourier transform 
which is difficult to evaluate directly. It is easier to start with the Fourier 
transform itself and work backwards using the inverse Fourier transform. 
Suppose we want to find the time-domain signal which has Fourier 
transform X (72) = 6 (2 — Qo). We can begin by using the inverse 
Fourier transform [link] 


Equation: 
1 * jt 
x(t) = — | 6(Q—M)eiaN 
2m J = 
_ 1 oa 
7 On « 


This result follows from the sifting property of the impulse function. By 
linearity, we can then write 
Equation: 


e)* 2, 25 (Q— MN) 


Cosine 


The cosine signal can be expressed in terms of complex exponentials using 
Euler's Identity 
Equation: 


cos (Qt) = (Ca + ev) 


bo] 


Applying linearity and the Fourier transform of complex exponentials to the 
right side of [link], we quickly get: 
Equation: 


cos (Mt) — 7d (MQ — Mo) + 7b (2+ MN) 


Real Exponential 


The real exponential function is given by x (t) = e~“u (t), where a > 0. 
To find its FT, we start with the definition 


Equation: 
_ . 
X(jQ) = / es 
0 
CO 
= / e (a+ IQ)t Gy 
0 
= FP _ 6 -(atsie| 
a+ 952 0 
—1 
= 0-1 
a+ 79 ( ) 
_ 1 
at 5 
therefore, 
Equation: 
1 


The Unit Step Function 


In a previous section, we looked at the unit step function, 


Equation: 
1, t>0 
t) = 7 
u(t) 1 t<0 


A closely related signal is the signum function, defined by 
Equation: 


1, t>0 
2 oe ee ee 


from which it follows that 
Equation: 


u(t) = (sgn(t) + 1)/2 
The signum function can be described as follows: 
Equation: 


sgn (t) =lim [e “u (t) — eu (—t)| 


a—0 


Since we already have the Fourier transform of the exponential signal, 
Equation: 


1 1 
SGN) =) — 
2) a0 | a +j2 a—jQ 
; —2522 
=|] 
a>0 q?2 -+- 22 
2 


Using [link] and linearity then leads to 


Equation: 


x(t) 


rect(t, T) 


u(t) + 15 (2) += 


X (522) 


sin(27/2) 
Qr/2 


1 
276(2 = No) 


10 (2 — 29) + 76 (2+ Mo) 


noe 
atjQ 


Some common Fourier transform pairs. 


When working problems involving finding the Fourier transform, it is often 
preferable to use a table of transform pairs rather than to recalculate the 
Fourier transform from scratch. Often, transform pairs in can be combined 
with known Fourier transform properties to find new Fourier transforms. 


Example 3.1 Find the Fourier transform of: y (t) = 2e°’u (—t). Clearly, 
we can write y(t) = x(—t) where a (t) = 2e~-°*u (t). Therefore, we can 


combine the known transform of x(t) from [link], namely, 
Equation: 


2 
5+ 502 


X (j2) = 


with the time reversal property: 
Equation: 


a (—t) & X(jQ)” 


to get the answer: 
Equation: 


Fourier Transform of Periodic Signals 


If the signal of interest is periodic with period 7’, then it has a Fourier 
Series: 
Equation: 


CO 
x(t) = S° Ces ort 


n=—OCo 


Using the linearity of the Fourier Transform, we have 
Equation: 


e.@) 


Ss GF {ei 


n=—OoO 


20 es Cn 0 (2 — nM) 


n=— OOo 


X(j22) 


where F'{ } corresponds to the Fourier Transform of the signal within the 
brackets. 


Filters 


Filters are devices which are commonly found in electronic gadgets. When 
you adjust the bass (low frequency) or treble (high frequency) settings on 
your MP3 player, you are adjusting the characteristics of a filter. A more 
technical name for a filter is a linear system . A filter is represented by a 
box having a single input (usually x(t)) and a single output (say, y(t)) as 
seen in [link]. 


— >| LI] 


Continuous-time 
filter. 


We can denote the operation the filter has on the input using the following 
notation: 
Equation: 


The types of filters we will consider in this book are linear and time- 
invariant. A filter is time-invariant if given that y(t) = L/a(t)], then 

y(t — rT) = L|x(t — 7). In other words, if the input to the filter is delayed 
by 7, then the output is also delayed by T. A filter is linear if given that 

yi (t) = Liay (t)] and yo (t) = L[ae (t)] then 

Equation: 


ay (t) + Bye (t) = Llaz, (t) + Br (t)| 


Equation [link] is often referred to as the superposition principle. We can 
use linearity and time invariance to derive the mathematical operation 
which the filter performs on the input, x(t). To do this we begin with the 
assumption that 


Equation: 


The signal h(t) is called the impulse response of the filter. From time 
invariance, we have 
Equation: 


h(t—7) = L[d(t — 7) 


Now we can use linearity to find the filter output when the input is 
x(7)d(t — 7), where x(7) is a constant 
Equation: 


x(r)h(t — 7) = L[a(r)d(t — 7)] 


We can extend the linearity property further by noting that 
Equation: 


where we can assume that the constants 7,, are ordered so that 

Ti, < Th, i < kand A, = T, — Tn—1. In [link], we are simply multiplying 
each 6(t — T,) by the constant x (7;,) An, so once again linearity should 
prevail. Now if we take the limit A, — 0, we obtain 

Equation: 


[ e@ne—nar=2| [eae nar 


CO [e,2) 


Using the sifting property of the unit impulse in the right side of [link] gives 
Equation: 


So it follows that the filter performs the following operation on the input, 


w(t): 


Equation: 


L{x(t)] 
- / 2 (r)h(t — r)dr 


CO 


y(t) 


The integral in [link] is called the convolution integral. A change of 
variables can be used to show that 
Equation: 


[. On eee ae a othe Dar 


CO —oo 


which means that the order in which two signals are convolved is 
unimportant. A short-hand notation for convolution is 
Equation: 


[- n\n adr=o ne) 


CO 


Properties of Convolution Integrals 
We list several important properties and their proofs. 


1. Commutative Property: 
Equation: 


z(t)*h(t) = h(t)*z(t) 
Lets start with 
Equation: 


x(t)*h(t) = : aes 


—CoO 


and make the substitution ~ = t — T. It follows that 
Equation: 


a(t}*n(e) =f 2(t—a)h(a)drar 


Ps 
—. 
oF 
aS 
* 
8 
—. 
oF 
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2. Associative Property: 
Equation: 


[a (t)*hi (t)|"ho (t) = a (t)*[hi (t)* he (6)] 


To prove this property we begin with an expression for the left-hand 
side of [link] 
Equation: 


/ Sone 
where we have expressed x (t)*h, (¢) as a convolution integral. 
Expanding the second convolution gives 

Equation: 


i; “ if GE "dr flea 


CO (ee) 


Reversing the order of integration gives 


Equation: 
[. x (T) if hi (y — T)he (t — Neh] oe 


CO CO 


Using the variable substitution ¢ = y — 7 and integrating over » in 
the inner integral gives the final result: 


Equation: 
[2] [_ m@rete—7- g)ar]ar 


CO CO 


where the inner integral is recognized as h, (t)*hg (t) evaluated at 
t = t — 7, which is required for the convolution with x(t). 

3. Distributive Property: 
Equation: 


a (t)*[ha (t) + he (€)] = @ (t)*ha (t) + @ (E)*he (t) 


This property is easily proven from the definition of the convolution 
integral. 

4. Time-Shift Property: If y(t) = x(t)*h(t) then 
x (t — to)*h (t) = y(t — to) Again, the proof is trivial. 


Evaluation of Convolution Integrals 


The key to evaluating a convolution integral such as 
Equation: 


CO 


OO / setae 


—oo 


is to realize that as far as the integral is concerned, the variable ¢ is a 
constant and the integral is over the variable 7. Therefore, for each t, we are 
finding the area of the product x(7)h(t — 7). Let's look at an example that 
illustrates how this works. 


Example 3.1 Find the convolution of x(t) = u(t) and h(t) = e*u(t). 
The convolution integral is given by 
Equation: 


[link] shows the graph of e~7u (rT), e~‘u (t), and their product. From the 
graph of the product, it is easy to see the the convolution integral becomes 


Equation: 
[evar= 1—e™, tee 
0 0, t<0 


Graphs of signals used in Example [link]. 
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S 
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Signals which can be expressed in functional form should be convolved as 
in the above example. Other signals may not have an easy functional 
representation but rather may be piece-wise linear. In order to convolve 
such signals, one must evaluate the convolution integral over different 
intervals on the t-axis so that each distinct interval corresponds to a 
different expression for x(t)*h(t). The following example illustrates this: 


Example 3.2 Suppose we attempt to convolve the unit step function 
z(t) = u(t) with the trapezoidal function 
Equation: 


f 02721 
Rijs dy lee? 
0, elsewhere 


From [link], it can be seen that on the interval 0 < t < 1, the product 

z(t — T)h(r) is an equilateral triangle with area ¢? /2. On the interval 

1 <t < 2, the area of x(t — r)h(r) is t — 1/2. This latter area results by 
adding the area of an equilateral triangle having a base of 1, and the area of 
a rectangle having a base of ¢ — 1 anda height of 1. For all values of t 
greater than 2, the convolution is 1.5 since x(t — r)h(r) = h(r) and h(r) 


is a trapezoid having an area of 1.5. Finally, for t < 0, the convolution is 
zero since x(t — T)h(r) = 0. 


Graphs of signals used in Example [link]. 
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Frequency Response 


Recall from [link] that the convolution integral 
Equation: 


has the Fourier Transform: 
Equation: 


¥(jQ) = H(jQ)X(jQ) 


where H(jQ) and X(j2) are the Fourier Transforms of h(t) and x(t), 
respectively. Solving for H(j) gives the frequency response: 
Equation: 


H (j2) = 


The frequency response, the Fourier Transform of the impulse response of a 
filter, is useful since it gives a highly descriptive representation of the 
properties of the filter. The frequency response can be considered to be the 
gain of the filter, expressed as a function of frequency. The magnitude of 
the frequency response evaluated at 2 = 2, |H(j{2o)| gives the factor the 
frequency component of x(t) at 2 = would be scaled by. The phase of 
the frequency response at 2 = 929, ZH(j{2) gives the phase shift the 
component of x(t) at 2 = Mp would undergo. This idea will be discussed 
in greater detail in [link]. A lowpass filter is a filter which only passes low 
frequencies, while attenuating or filtering out higher frequencies. A 
highpass filter would do just the opposite, it would filter out low 
frequencies and allow high frequencies to pass. [link] shows examples of 
these various filter types. 


re) 
(a) 


0 
JHGQ)| 
(b) 
0 Q 
|HGQ)| 
(c) 
7 Q 
Different filter types: 


(a) lowpass, (b) 
bandpass, (c) highpass. 


The Sinusoidal Steady State Response 


It is useful to see what the effect of the filter is on a sinusoidal signal, say 
x (t) =cos (Mot). If y(t) is the output of the filter, then we can write 
Equation: 


woe / Ee aide 


CO 


Using the Euler formula for cos ({2ot), right hand side of [link] can be 
written as: 
Equation: 


1 / ~ (citoute—r) 7 e-natt-7)) nae 
2 Jo 


This integral can be split into two separate integrals, and written as: 
Equation: 


GNot oe) jot CO 
: ; e I%Th (r)dr + : A / el Th (1) dr 


2 
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The first of the two integrals can be recognizes as the Fourier Transform of 
the impulse response evaluated at {2 = §29. The second integral is just the 
complex conjugate of the first integral. Therefore [link] can be written as: 
Equation: 


eI Mot . e IMot 
H (jo) + 


H™ (jo) 


Since the second term in [link] is the complex conjugate of the first term, 
we can express [link] as: 
Equation: 


Ref{el"H (jo) } 


or expressing H(j2o) in terms of polar coordinates: 
Equation: 


Ref e3 |H (j%p)\e'#0%) | = Ref [A (jp) [estat <HU500) | 


Therefore, we find that the filter output is given by 
Equation: 


y(t) = | (9%) |cos (Qot + ZH(GN)) 


This is called the sinusoidal steady state response. It tells us that when the 
input to a linear, time-invariant filter is a cosine, the filter output is a cosine 
whose amplitude has been scaled by |H(j§29)| and that has been phase 
shifted by ZH (729). The same result applies to an input that is an 
arbitrarily phase shifted cosine (e.g. a sine wave). 


Example 3.1 Find the output of a filter whose impulse response is 

h (t) = eu (t) and whose input is given by x(t) =cos (2t). It can be 
readily seen that the frequency response of the filter is 

Equation: 


1 


se ae rary) 


and therefore |H(72)| = 0.1857 and 7H(j2) = —0. 3805. Therefore, 
using [link]: 
Equation: 


y(t) = 0.1857 cos (2t — 0. 3805) 


Parallel and Cascaded Filters 


In some applications, such as graphic equalizers, it is useful to place filters 
in parallel as shown in [link]. Can the parallel combination of filters be 
characterized by a single equivalent filter he, (t)? The answer is yes and 
results by noting that 

Equation: 


y(t) = Sa (t)*hi(t) 


= [-2t-n omer 


Therefore, the last equation in [link] shows that 
Equation: 


S| tt 
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Parallel filter structure. We wish to find 
an equivalent filter with impulse 
response he, (t). 


The equivalent transfer function for the parallel filter structure is given by 
Equation: 


eq (j22) — A; (j22) 


Me 


I 
=" 


t 


Next we wish to find an equivalent filter for the cascaded structure shown 
in [Link]. 


t t 
x(t) hi yi(t) y(t) 
AO >| hat 


Cascaded filter structure. We wish to 
find an equivalent filter with impulse 
response Neg (t). 


This can be done by finding an expression for the intermediate value yj (t): 
Equation: 


= / seine: 


[o,@) 


The output of the cascaded structure is given by 
Equation: 


y= f  He=AGe 


o,8) 


substituting [link] into [link] gives 
Equation: 


yit)= [| f ete—9— 2m (r)ar| ha Nan 


Reversing the order of integration and rearranging slightly gives 
Equation: 


y= ff 2(¢—4=n)ha (r)ha (hey 


Now let € = y + 7, solving for 7 gives 7 = € — y and dé = dr. 
Substituting these quantities into [link] leads to 
Equation: 


v= ae-ol f ne =a 


Notice that we can factor x(t — €) from the inner integral since x(t — €) 
does not depend on +. The integral in the brackets is recognized as 

hy (t)*hg (t) evaluated at €. Therefore for the cascaded system, the 
equivalent impulse response is given by 

Equation: 


het) = f ” ha (¢— Nha DAY 


oo 


This can be generalized to any number of cascaded filters giving 
Equation: 


Neg (t) = ha (t)*hi (t)*: - -*Ay (é) 


First Order Filters 


A first-order lowpass filter has the frequency response 
Equation: 


1 


Hyp (G2) = —— 
1+ j5- 


The frequency at which the frequency response magnitude has dropped to 
1/ V2 is called the corner frequency[footnote]. The frequency response 
magnitude and phase are plotted in [link]. It is common to express the 
frequency response magnitude in units of decibels (dB) using the formula 
This term most likely originates from it's role in Bode plots, a shortcut 
method for sketching the graph of a frequency response. 

Equation: 


20logio| HA (j2)| 


At the corner frequency for a first order lowpass filter, the frequency 
response magnitude is 1/ 4/2 or roughly -3 dB. From [link], it can easily be 
seen that the impulse response for the first-order lowpass filter is given by 
Equation: 


hyp (t) = Ne Pty (t) 


A first-order highpass filter is given by 
Equation: 


‘ 2. 
Hyp (jQ) = —— 
1-7 


Notice that 
Equation: 


1 


Hyp (jf) = 1 — ——— 
1+ je 


This makes sense since a highpass filter can be constructed by taking the 
filter input x(t) and subtracting from it a lowpass filtered version of x(t). 
The impulse response of the first-order highpass filter therefore becomes: 
Equation: 


hip (t) = 6(t) — Qe u(t) 


IHGQ)! (dB) 


ZHGQ) (rad) 
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Frequency response magnitude and phase for a first-order 
lowpass filter ((2, = 1 rad/sec). 


IHGQ)! (dB) 


ZHGQ) (rad) 
= 
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Frequency response magnitude and phase for a first-order 
highpass filter (2, = 1 rad/sec). 


First order filters can be easily implemented using linear circuit elements 
like resistors, capacitors, and inductors. [link] shows a first order filter 
based on a resistor and a capacitor. Since the impedance for a resistor and 
capacitor are R and 1/j2C;, respectively, voltage division leads to a 
frequency response of 

Equation: 


1 


H 2) = ———__—_ 
a (J ) 14+ 72R\C, 


Therefore the corner frequency for this filter is (2, = roe Similarly, a 
first-order highpass filter can be implemented using a resistor and capacitor 


as shown in [link]. This filter has a frequency response of 
Equation: 


j QR» Cy 


H 2 => — 
up(J ) 1 + jQR C3 


The corner frequency for the highpass filter is seen to be 2, = 


Ry 


+ 
v;(t) C l Vo(t) 


Circuit implementation of a 
first-order lowpass filter 
having Q, = 1/RiCi. 


C, 


+f 
v;(t) Ry V(t) 


Circuit implementation 
of a first-order highpass 
filter having 
W, ol / RoC. 
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Now one might be tempted to apply the results of [link] to build a bandpass 
filter by cascading the lowhpass and highpass circuits in Figures [link] and 
[link], respectively. Theory would predict that the equivalent frequency 
response of this circuit is given by 

Equation: 


Fleg (92) = Hyp (G2) Hyp (52) 


Unfortunately, this is not possible since the circuit elements in the lowpass 
and highpass filters interact with one another and therefore affect the 
overall behavior of the circuit. This interaction between the two circuits is 
called loading will be studied in greater detail in the exercises. To get 
theoretical behavior, it is necessary to use a voltage follower circuit, 
between the lowpass filter from the highpass circuits. The voltage follower 
circuit is usually an active circuit (requires external power supply) that has 
very high input impedance and very low output impedance. This eliminates 
any loading effects which would normally occur between the lowpass and 
highpass filter circuits. 


Parseval's Theorem for the Fourier Transform 


In Chapter 2, we looked at a version of Parseval's theorem for the Fourier 
series. Here, we will look at a similar version of this theorem for the Fourier 
transform. Recall that the energy of a signal is given by 

Equation: 


If the energy is finite then x(t) is an energy signal, as described in Chapter 
1. Suppose x(t) is an energy signal, then the autocorrelation function is 
defined as 

Equation: 


Ty (¢) = a (t)*e (—t) 


It can be shown that r, (t) is an even function of t and that r, (0) = e,(see 
Exercises). The Fourier transform of r, (t) is given by 


X (jQ)X(jQ) = |X(jQ)|°. If follows that 


Equation: 
SAS Ae ” 1xX(GQ) Pela 
‘ 27 Jo t=0 
1 - eo 
= — X(j2Q)|"dQ 
a | MIP 


Which is Parseval's theorem for the Fourier transform. 


The Fourier Transform: Excercises 


Exercises 
1. Find the Fourier transform of the following signals. Sketch the graph 


of the Fourier transform if it is real, otherwise, sketch the magnitude 
and phase fo the Fourier transform: 


x ( 
x(t) = 4e°7y (—t) 
x ( 


7a nwoan oo ow 
8 
—_ —™ 
Cs 
Ne: 
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a(t) = rect(t, 1) 
et) Ae we 
x(t) = cos(10zt) 
a(t) =6 
eee i jt] < 0.5 
OF jé| > 0.5 


2. Find the convolution of the following pairs of signals: 


x(t) h(t) 


rd A 
t t 
a. 0 1 0 0.5 


x(t) h(t) 


rd 1. 
t 
b. ° 0 


3. Find the output of the filter whose transfer function is 
Equation: 


10. 


21 


I) eG 


and whose input is z(t) = u(t). Hint, find the impulse response h(t) 
corresponding to H(792) and convolve it with the input. 


. Show that if v (t) = D[u(t)], then 


Equation: 


i. "gD eS i i a] 


Hint: Integrate both sides of v (t) = L/u(t)]. Then express the right 
hand integral as the limit of a sum (as in a calculus textbook). Then by 
linearity, you can exchange the sum and the L]-]. 


. Find an expression for the convolution of x(t) = u(t) and 


h(t) = sin(8t)u(t) 


. Find an expression for the convolution of x(t) = rect(¢ — 0.5, 1) and 


h(t) =e “a (ft). 


. Find the Fourier transform of the periodic signal in problem 2, Chapter 


2. 


. Consider a filter having the impulse response 


Equation: 
h(t) =e “u(t) 


Sketch the frequency response (both magnitude and phase) of the filter 
and find the output of the filter when the input is x(t) =cos (10¢). 


. Repeat the previous problem for the impulse response given by 


Equation: 


lL. OSA 
Oe. ans 
(1) a otherwise 


Suppose that two filters having impulse responses h (t) and hg (t) are 
cascaded (i.e. connected in series). Find the impulse response of the 
equivalent filter assuming hy (t) = 10e~'°'w (t) and 

hy (t) = 5e- u(t). 


11. Design a first-order lowpass filter having a corner frequency of 100 
Hz. Use a 100K? resistor. Plot both the magnitude and phase of the 
filter's frequency response. 

12. Design a first-order highpass filter having a corner frequency of 1000 
Hz. Use a 0. O1yF capacitor. Plot both the magnitude and phase of the 
filter's frequency response. 

13. The following problems are associated with the circuits in [link]: 


0.01uF 
(a) uke (b) 
: o 
v;(t) IuF == v(t)  vi(t) SkQ> V(t) 
aA 0.01uF 
(c) . 
v;(t) IMF ska> v,(t) 
Problem [link]. 


a. Find the frequency response of the circuit in [link](a), and sketch 
its magnitude and phase. 

b. Find the frequency response of the circuit in [link](b) and sketch 
its magnitude and phase. 

c. Find the frequency response of the filter in [link](c), sketch its 
magnitude and phase and show that it is not the product of the 
frequency responses for problems [link] and [link]. 


The Laplace Transform: Introduction 
Introduces the double-sided and single-sided Laplace transforms. 


The Double-Sided Laplace Transform 


You may have noticed that we avoided taking the Fourier transform of a 
number of signals. For example, we did not try to compute the Fourier 
Transform of the unit step function, z(t) = u(t), or the ramp function 

x(t) = tu(t). The reason for this is that these signals do not have a Fourier 
transform that converges to a finite value for all (2. Recall that in order to 
deal with periodic signals such as cos (Qt) we had to settle for Fourier 
transforms having impulse functions in them. The Laplace transform gives 
us a mechanism for dealing with signals that do not have finite-valued 
Fourier Transforms. 


The double-sided Laplace Transform of x(t) is defined as follows: 
Equation: 


Xie) = / ” 2 (teat 


OO 


where we define 
Equation: 


s=ot+ gj 


We observe that the Laplace transform is a generalization of the Fourier 
transform since 
Equation: 


X (G2) = X(8)|.—-j0 


Therefore, we can write 
Equation: 


X (s) = F{e “x (t)} 


The inverse Laplace transform can be derived using this idea. Applying the 
inverse Fourier transform to F{e~°‘z (t) } gives 


Equation: 


e “x (t) : [- Fle “a (t) bed 


Ot dc 


Using [link] leads to 
Equation: 


dil = — 


Substituting s for 2 and solving for x(t) in [link] gives 
Equation: 
1 o+joo 
x (t) X (s)e"ds 


(On a 


Equation [link] is called the inverse Laplace transform. The integration is 
along a straight line in the complex s-plane corresponding to a fixed value 
of o. This is illustrated in [link]. 


The inverse Laplace transform 
integrates along a line having a 
constant o in the complex s— 
plane. 


It is important that this line exist in a region of the s-plane that corresponds 
to the region of convergence for the Laplace transform. The region of 
convergence is defined as that region in the s-plane for which 

Equation: 


ii |x (t)e **|dt < 00 


CO 


Note that since s = 0 + j£2 this is equivalent to 
Equation: 


/ |x (t)e “|dt < 00 


e.@) 


The Single-Sided Laplace Transform 


We define the single-sided Laplace transform as 
Equation: 


R= i ” 2 (t)e~*dt 


where the lower limit of integration tacitly includes the point t = 0. That is, 
0” represents a point just to the negative side of ¢ = 0. This allows for the 
integral to take into account signal features that occur at t = 0 such asa 
step or impulse function. The single-sided Laplace transform is motivated 
by the fact that most signals are turned on at some point. If the region of 
convergence is not specified then different signals can yield identical 
double-sided Laplace transforms, as the following example illustrates. 


Example 3.1 Consider the signal 
Equation: 


z(t) =e “u(t) 


The double-sided Laplace transform is given by 
Equation: 


X1(s) = / e “e “dt 
0 


= / o (stat ge 
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I] 0 
=— ep (tj 2+a)t 


s+a 0 


— =n a ee = 1) 
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the magnitude of e~ ‘7+ is zero only if ¢ > —a which establishes the 


region of convergence. Therefore we have 
Equation: 


e“u(t)o 


,o > -a 
s+a 


Now consider the Laplace transform of the signal 


Equation: 
z(t) = —e u(t) 
We have 
Equation: 
0 
Yale) = —f eteta 
CO 
0 
— -f « —(s+ta t dt 
CO 
— 1. -oriarays 
So os 
= 1 (1 = e(7ta)joo ee) 
S+aQ 


Here, the quantity e(°+@) is zero only if ¢ < —a so we have 
Equation: 


—e “u(-t)o so < -a 


s+a 


which is identical to X (s) except for the region of convergence. 


To avoid scenarios where two different signals have the same double-sided 
Laplace transform, we restrict our signals to those which are assumed to be 
zero for t < 0. Such signals are called causal signals. 


The region of convergence for the single-sided Laplace transform is a 
region in the s-plane satisfying o > Omjin as shown in [link]. 


jQ 


The region of convergence of the 
single-sided Laplace transform in 
the complex s—plane. 


To see this, we observe that if 


Equation: 
[e,@) 
| |x (t)e °m"" dt < 
ae 


then it must be the case that 
Equation: 


| |x (t)e *'|dt < 00 


for > Omin, since e “ decreases faster than e~ 7n*. Finally, the inverse 


single-sided Laplace transform is the same as the inverse double-sided 
Laplace transform (see [link]), since a single sided Laplace transform can 
be interpreted as the double-sided Laplace transform of a signal satisfying 
x(t) = 0,t < 0. From here on, we will work exclusively with the single- 
sided Laplace transform. Unless we need to specifically differentiate 
between the single or double-sided transforms, we will refer to the single- 
sided Laplace transform as simply the “Laplace transform". 


Properties of the Laplace Transform 
Derives basic properties of the Laplace transform. 


Properties of the Laplace Transform 


The properties associated with the Laplace transform are similar to those of the Fourier 
transform. First, let's set define some notation, we will use the notation £{} to denote the 
Laplace transform operation. Therefore we can write X (s) = £{x(t)} and 

ax (t) = £~'{X(s)} for the forward and inverse Laplace transforms, respectively. We can also 
use the transform pair notation used earlier: 

Equation: 


z(t) + X(s) 


With this notation defined, lets now look at some properties. 


Linearity 


Given that x, (t) <> Xj (s) and x2 (t) + X2 (s) then for any constants a and (, we have 
Equation: 


ax, (t) + Bro (t) 6 aX, (s) + BX2 (s) 


The linearity property follows easily using the definition of the Laplace transform. 


Time Delay 


The reason we call this the time delay property rather than the time shift property is that the time 
shift must be positive, i.e. if 7 > 0, then x(t — 7) corresponds to a delay. If r < 0 then we 
would not be able to use the single-sided Laplace transform because we would have a lower 
integration limit of 7, which is less than zero. To derive the property, lets evaluate the Laplace 
transform of the time-delayed signal 

Equation: 


eee / "eae ea 


Letting y = t — 7 leads tot = y +7 and dt = dy. Substituting these quantities into [link] gives 
Equation: 


L{x(t —7)} 


I| 
(ieee 
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ee 
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=) 


0 (oe) 
ae | x (ye dy + a xz (ye "dy 
- 0 


where we note that the first integral in the last line is zero since x(t) = 0,t < 0. Therefore the 
time delay property is given by 
Equation: 


L{a(t-—7T)} =e "X(s) 


s-Shift 

This property is the Laplace transform corresponds to the frequency shift property of the Fourier 
transform. In fact, the derivation of the s-shift property is virtually identical to that of the 
frequency shift property. 

Equation: 


Se Va yt = i e “2s (the “dt 
0 

= a (t)e@*9)'dt 
0 


= i x (te! 


= X(s+a) 


Q 


+o4jQ)t dt 


The s-shift property also alters the region of convergence of the Laplace transform. If the region 
of convergence for X(s) is o > Omin, then the region of convergence for Bieree (t)} is 


Multiplication by ¢ 


Let's begin by taking the derivative of the Laplace transform: 
Equation: 


dX(s) af © int 
= — the “dt 
ds ds Jo A 


me d 
— je “a 
i x (t)—~e 


s 


_ / ta (t)e~ ‘dt 
0 


So we can write 
Equation: 


dX(s) 
ds 


S{ta(t)} = — 


This idea can be extended to multiplication by t”. Letting y(t) = ta(t), if follows that 
Equation: 


dY(s) 
ty(t) ~ — 
y(t) as 
d?X 
4 EX) 
ds? 
Proceeding in this manner, we find that 
Equation: 
n a” X (s) 


Time Scaling 


The time scaling property for the Laplace transform is similar to that of the Fourier transform: 
Equation: 


L{x(at)} = [PO a(atye a 


1s : 
= —|[ ax(qyJe ody 
a Jo 
Zea x(+) 
a a 


‘ : aches d 
where in the second equality, we made the substitution t = 2 and dt = a 


Convolution 


The derivation of the convolution property for the Laplace transform is virtually identical to that 
of the Fourier transform. We begin with 
Equation: 


af Ce dr} 7 [- se nea ar 


CO CO 


Applying the time-delay property of the Laplace transform gives 
Equation: 


iP “ORG ie > He / ” ete 
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If h(t) is the the impulse response of a linear time-invariant system, then we call H(s) the 
system function of the system. The frequency response results by setting s = j2 in H(s). The 
system function provides us with a very powerful means of determining the output of a linear 
time-invariant filter given the input signal. It will also enable us to determine a means of 
establishing the stability[ footnote] of a linear-time invariant filter, something which was not 
possible with the frequency response. 

We will discuss stability shortly 


Differentiation 


The Laplace transform of the derivative of a signal will be used widely. Consider 


Equation: 
eS a / ” 2! (t)e-*dt 
dt 2 


this can be integrated by parts: 
Equation: 


which gives 
Equation: 
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therefore we have, 
Equation: 


Higher Order Derivatives 


The previous derivation can be extended to higher order derivatives. Consider 
Equation: 


ga EO ei) 
dt 

it follows that 
Equation: 

dy(t) 

— Y (s) —y(0 

a (s)—y(0") 
which leads to 
Equation: 
- dzx(0-) 


This process can be iterated to get the Laplace transform of arbitrary higher order derivatives, 
giving 
Equation: 


where it should be understood that 


Equation: 
d™z(0-) d™x(t) 
———__—_ = »m=]1, 
dt™ Oe is 
Integration 
Let 
Equation: 


it follows that 


Equation: 
dg(t) 
SNEED a t 
a Te) 
and g(0) = 0. Moreover, we have 
Equation: 
dg(t) 
ae) 0 | 
= sG(s)—g(0-) 
= sG(s) 
therefore 
Equation: 
X(s 
G(s) = ) 
but since 
Equation: 


we have 
Equation: 


Now suppose z(t) has a non-zero integral over negative values of t. We have 


Equation: 
i. eae [. ere [ ear 


The quantity f oe x (r)dr is a constant for positive values of t, and can be expressed as 
Equation: 


it follows that 
Equation: 


where we have used the fact that u(t)  -. 


The Initial Value Theorem 


The initial value theorem makes it possible to determine x(t) at t = 0* from X(s). From the 
derivative property of the Laplace transform, we can write 


Equation: 
ef ae) \ = sX(s) —2(0°) 


dt 


Taking the limit s — oo 
Equation: 


lim * GEN): cdg = lim [sX (s) _ «(0-)| 


800 Jy_ dt 800 


Os. det), -. ; : 
[ Jim a © ‘dt =lim [sX (s) — «(0-)| 


There are two cases, the first is when x(t) is continuous at t = 0. In this case it is clear that 
Gt) est — 0 as s — o, so [link] can be written as 
Equation: 


0 =lim [sX (s) — 2(0-)| 


Since x(t) is continuous at t = 0, x(0~) = (0°), the Initial Value Theorem follows, 
Equation: 


z(0T) =lim sX(s) 


s—0o 


The second case is when x(t) is discontinuous at t = 0. In this case, we use the fact that 
Equation: 
dx(t) 
dt 


= [2(0") — 2(0°)]6@) 


t=0 


For example, if we integrate the right-hand side of [link] with z(0~) = 0 and x(0*) = 1, we 
get the unit step function, u(t). Proceeding as before, we have 
Equation: 


The left-hand side of [link] can be written as 
Equation: 


07 lo) 


lim [ [2(0*) —2(0-)] eats tim f° PO ese 


s—0o ot s—0o ot 


From the sifting property of the unit impulse, the first term in [link] is 
Equation: 


[z(0") —2(0°)| 


while the second term is zero since in the limit, the real part of s goes to infinity. Substituting 
these results into the left-hand side of [link] again leads to the initial value theorem, in [link]. 


The Final Value Theorem 


The Final Value Theorem allows us to determine 
Equation: 


lim z (t) 


t-0o 


from X(s). Taking the limit as s approaches zero in the derivative property gives 
Equation: 
CO 


: dz(t) _, : 7 
lim ae ‘dt =lim [sX (s) — 2(0-)| 


s—0 0- 


The left-hand-side of [link] can be written as 
Equation: 


[ lim OO) sat = a aly) = 2 (oo) — (07) 


_ s>0 - dt 


Substituting this result back into [link] leads to the Final Value Theorem 
Equation: 


x (oo) =lim sX (s) 


s—0 


which is only valid as long as the limit z(0o) exists. 


Property y(t) Y(s) 
Linearity ax, (t) + Bre (t) aX 1(s) + BX2(s) 
Time Delay z(t — T) X (sje 7 


s-Shift xz (t)e~* X(s +a)) 


Multiplication 
by t 


Multiplication 
by t” 


Convolution 


Differentiation 


Integration 
Initial Value 
Theorem 


Final Value 
Theorem 


tz(t) _ dX(s) 
ds 
t"x (t) (—1)" ees) 
x(t)*h(t) X(s)H(s) 
dx(t) = s(o- 
- sX (s) (0-) 

dx (t) ‘ - dzx(07) 
1 s°X (s) — sx(0-) — a 
da (t) s"X(s)—s" ‘x(0-) — : 8 
= X(s) OS De 


Laplace Transform properties. 


Laplace Transforms of Common Signals 
This module derives the Laplace transform of well-known signals. 


Laplace Transforms of Some Common Signals 


We'll next build up a collection of Laplace transform pairs which we will 
include in a table. It's important to keep in mind that once the transform pair 
has been derived, the focus should be on utilizing the transform pair found 
in the table rather than in recalculating the transform. 


Exponential Signal 


Consider the Laplace transform of x (t) = e“u (t): 


Equation: 
CO 
L{e“u(t)} = / e “edt 
0 
ee) 
= / oe (atsyt ay 
0 
— at o-(atorjat|™ 
a+s 0 
a Ae ee 


where o > —a defines the region of convergence. Notice also that if a < 0 
, X(s) still exists provided o > —a. Therefore, 
Equation: 


Unit Step Function 


Recall that the Fourier transform of u(t) does not converge for this signal 
(it contains a 6(2) term, since u(t) is not absolutely integrable). The 
Laplace transform on the other hand, converges. In fact, since u(t) is a 
special case of the exponential function with a = 0, 

Equation: 


The region of convergence is ao > 0, the right-half plane. 


Ramp Signal 


This signal is given by x(t) = tu(t) and also does not have a finite Fourier 
transform. The Laplace transform can be easily found using the 
multiplication by ¢ property, 

Equation: 


The region of convergence is once again the right-half plane, o > 0. 


Cosine Signal 


Even though we computed the Fourier transform of the cosine signal, 

x (t) =cos (Mot), the Fourier transform technically does not converge for 
this signal. That is why X (j2) involves impulse functions. The Laplace 
transform produces quite a different result. First we use the fact that 
Equation: 


eF@oty (t) + e Ft (t) 


cos (Qot)u (t) = 5 


Since each of the two terms is an exponential function we have 
Equation: 


L{cos (Npt)u(t)} = 
s2+2? 


Here, the region of convergence corresponds to the right-half plane, 0 > 0. 


More Transform Pairs 


We can use the existing transform pairs along with the properties of the 
Laplace transform to derive many new transform pairs. Consider the 
exponentially weighted cosine signal. This signal is given by 
Equation: 


a(t) =e “ cos (Qot)u(t) 


We can use the s—shift property of the Laplace transform [link] along with 
the Laplace transform of the cosine signal [link] to get 
Equation: 


e * cos (Qt)u(t) o = 
(sta)°+ 22 


Another common signal is 
Equation: 


az (t) =te “u(t) 


We use the Laplace transform of the exponential signal !"Exponential 
Signal" and the ¢ multiplication property to get 
Equation: 


1 
te “u(t) o —_— 
(s +a) 


Extending this idea one step further, we have 
Equation: 


az (t) = t?e-“u (t) 


Here, multiplication by t? applies, giving 
Equation: 

2 
te “ul(t)o ae 
(s+ a) 


Example 3.1 Consider the signal x (t) = te -*“u (t). Therefore, we get 


Equation: 

1 
te u(t) o = 
(s + 2) 


Example 3.2 Consider the signal x (t) = e~*' cos (5t)u (t). As seen in 
[link], 
Equation: 


s+2 


e * cos (5t)u ce ae ENS 


cos (Qot)u (t) 


sin (Qot)u (t) 


e~™ cos (Qot)u (t) 


e ™ sin (Mot)u (t) 


teu (t) 


teu (t) 


t"e “uy (t) 


Some common Laplace Transform pairs. 


Finding the Inverse Laplace Transform 


Finding the Inverse Laplace Transform 


Using Transform Tables 
The inverse Laplace transform, given by 
Equation: 
1 o+joo 
z(t) = — X (s)e"ds 
2mj o— JOO 


can be found by directly evaluating the above integral. However since this 
requires a background in the theory of complex variables, which is beyond 
the scope of this book, we will not be directly evaluating the inverse 
Laplace transform. Instead, we will utilize the Laplace transform pairs and 
properties . Consider the following examples: 


Example 3.1 Find the inverse Laplace transform of 
Equation: 


By looking at the table of Laplace transform properties we find that 
multiplication by e~!°* corresponds to a time delay of 10 sec. Then from 
the table of Laplace transform pairs , we see that 

Equation: 


1 
s+5 


corresponds to the Laplace transform of the exponential signal eu (t). 
Therefore we must have 


Equation: 


a (t) = e510) y (¢ — 10) 


Example 3.2 Find the inverse Laplace transform of 


Equation: 
1 
X (s) = —— 
(s + 2) 


First we note that from the table of Laplace transform pairs , the Laplace 
transform of tu(t) is 
Equation: 


Then using the s-shift property in the table of Laplace transform properties 
gives 
Equation: 


az (t) = te **u (t) 


Also, the same answer may be arrived at by combining the Laplace 
transform of e~?¢u (¢) with the multiplication by t property. 


Partial Fraction Expansions 


Partial fraction expansions are useful when we can express the Laplace 
transform in the form of a rational function, 
Equation: 


bys? + bg-18t 1 +.---+6;5+ bo 
ApsP + Gp 1s? 1 +-+--+a18 +a 
B(s) 
A(s) 


A rational function is a ratio of two polynomials. The numerator 
polynomial B(s) has order q, i.e., the largest power of s in this polynomial 
is q, while the denominator polynomial has order p. The partial fraction 
expansion also requires that the Laplace transform be a proper rational 
function, which means that q < p. Since B(s) and A(s) can be factored, 
we can write 

Equation: 


(s — 81) (s — B2)-++ (8 — Ba) 


ease eee a) 


The 8;,7 = 1, 2,...,q are the roots of B(s), and are called the zeros of 
X(s). The roots of A(s), are a;,7 = 1,...,p and are called the poles of 
X(s). If we evaluate X(s) at one of the zeros we get 

X(8;) = 0,i = 1,...,q. Similarly, evaluating X(s) at a pole 

gives| footnote] X(a;) = +oo,i = 1,..., p. The partial fraction expansion 
of a Laplace transform will usually involve relatively simple terms whose 
inverse Laplace transforms can be easily determined from a table of 
Laplace transforms. We must consider several different cases which depend 
on whether the poles are distinct. 

The actual sign would need to be evaluated at some value of s that is 
sufficiently close to the pole. 


Distinct Poles: 


When all of the poles are distinct (i.e. a; A a;,7 A 7) then we can use the 
following partial fraction expansion: 
Equation: 


A A A 
ee ee ee 
$—Q, $—Q2 S—Qp 


X (s)= 


The coefficients, A;,2 = 1,...,p can then be found using the following 
formula 
Equation: 


A; = X(s) (8 — ai) |, 9,54 = 1-5 


Equation [link] is easily derived by clearing fractions in [link]. The inverse 
Fourier transform of X(s) can then be easily found since each of the terms 
in the right-hand side of [link] is the Laplace transform of an exponential 
signal. This method is called the cover up method. 


Example 3.3 Find the inverse Laplace transform of 


Equation: 
2s — 10 
X(s) = —————_ 
(s) s?+38s4+2 
a 2s —10 
7 (s + 1)(s + 2) 
Since the poles are aj = —1 and a2 = —2 are distinct, we have the 
expansion 
Equation: 
A A 
X(s) = 1 2 
s+1l s+2 


Using [link] then gives: 
Equation: 


_ 28s —10 
s+2 s=—1 
= —12 
and 
Equation: 


A, = X(s)(s+ 2)|,—-2 


_ 2s — 10 
7 s+1 s=—2 
= 14 
Therefore, we get: 
Equation: 
—12 14 
(s) = 
s+1 s+2 


The inverse Laplace transform of X(s) can be found by looking up the 
inverse transform of each of the terms in the right-hand side of [link] giving 
Equation: 


az (t) = —12e tu (t) + 14e~**u (t) 


Repeated Poles: 


Let's consider the case when each pole is repeated, 
Equation: 


as B(s) 
ae) (s — ay)"'(s — ag)” +--+ (s —a,)”" 


where p; + po +--:+p, = p. In this case the partial fraction expansion 
goes like this: 


Equation: 
Ai Ai? Ay 
X(s) = ——— + — 4+ + 
S$—Q1  (s—a}) (s — a4) 
Ao Ag» Ad», 
Sa as + Savss + ee 
S—@2  (s—ag) (s — a2) 
+ ae 
Ar A,2 Ary: 
+ 4 4 + 
5— Or (s—ar,) (s — ay) 


We'll look at two methods. In the first method, the coefficients can be found 
using the following formula 
Equation: 


where 7 = 1, 2,...,7, k = 0,1,...,p; — 1 and 
Equation: 


Xj (s) = X(s)(s — ai)" 
Note that the computation of A; ,, does not require any differentiation, 
since k = 0. 


Example 3.4 Find the inverse Laplace transform of 
Equation: 


x@= = 1 


(s +2)’ 
Here we have a single repeated pole at s = —2. The expansion is therefore 
given by 
Equation: 
A A 
xG)= 11 4 13 


s+2  (s+2)° 


Using [link], we begin with k = 0 which corresponds to 


Equation: 
Ai2 =X (s)(s 7 2)°| 
= | ee 
= —3 
Next, we set k = 1 in [link] 
Equation: 
Aya = [X(s)(s +2)" 
= — S)(s 
oe ds ba 
d 
= —|s—1 
ds ee fess 
=i 


The partial fraction expansion is then given by 
Equation: 


Therefore, 
Equation: 


az (t) =e **u (t) — 3te~*u (t) 


In the second method, the coefficients A;,,,2 = 1,...,r can be found via 
the cover up method. The remaining coefficients, 

Ak.p;,t = 1,...,7, & = 1, ...,p; — 1 can be found by substituting values of 
s that are not equal to one of the poles in [link]. This leads to a system of 
linear equations which can be used to solve for the remaining coefficients. 
This method is generally preferable if the order of each repeated pole as 
well as the number of poles is sufficiently small so that the number of 
unknown coefficients is at most two for hand calculations. 


Example 3.5 Find the inverse Laplace transform of: 
Equation: 

8 
(s+1)° 
Ail ‘ Aj2 7 Ai3 ; 
stl (s+1) (s +1) 


X(8) = 


Using the cover-up method we can find Aj 3 as follows 
Equation: 


Ai3 = s|,-.=—1 


So we are left with 
Equation: 


sae 
‘Cosy 
Ait n Ai2 1 


s+l  (s+1)? (s+1)° 


Setting s = 0 in [link] leads to 


Equation: 

Aji t+Ai2=1 
and setting s = —2 in [link] gives 
Equation: 

—Aji + Ai2=1 


These choices of s were used to simplify the linear equations to the greatest 
extent possible. The solution to [link] and [link] is easily found to be 

Aj = 0 and Aj 2 = 1. The partial fraction expansion is given by 
Equation: 


xs 


(s+1)* (s+1) 


Using the corresponding Laplace transform pairs leads to 
Equation: 


x (t) = te-*u (£) — ste tu (t) 


Distinct and Repeated Poles: 


If a Laplace transform contains both distinct and repeated poles, then we 
would combine the expansions in [link] and [link]. Perhaps the easiest way 
to indicate this is by way of an example: 


Example 3.6 Find the inverse Laplace transform of 
Equation: 
s+2 
(s +1) (s+3)(s+5)° 
A A A A 
= 1, 42°, 431 3,2 
S+1 8+3 8s+5 (545) 


X(s) = 


The coefficients corresponding to the distinct poles can be found using 
[link]: 
Equation: 


Equation: 
A, = X(s)(s+3)|,—-3 
s+2 
(s+1)(s+5)? 
i 


s=—3 
8 
The coefficient Az corresponding to the double pole at s = —5 can be 


found using [link] with k = 0: 
Equation: 


A32 =X(s)(s+5)"|_ 
_ s+2 
~ (s+1)(s +3) 
—3 


8 


s=—5 


The remaining coefficient, A3 1 can be found using [link] with k = 1: 
Equation: 
d 


A3i; = —~X(s)(s +5)" 
ds 


s=_— 


-<| s+2 | 

ds | (s+1)(s+3) }|,__5 

(s? + 4s + 3) — (s + 2) (28 +4) 
(s? + 4s +3)? 


—95 
32 


Alternately, A3 1 can be computed by substituting the values obtained for 
Aj, A2 and A3 2 back into [link] and then substituting an arbitrary value for 
s that does not equal one of the poles as indicated earlier, like s = 0. This 
leads to a simple equation whose only unknown is A3 1. The partial fraction 
of X(s) is then given by: 


Equation: 
2 
X(s) = ——— 
(s +1)(s+3)(s+5) 
1 1 5 3 
32 8 32 8 


Applying the inverse Laplace transform to each of the individual terms in 
[link] and using linearity gives: 
Equation: 


2()= ae (t) ++ —e u(t) — —e u(t) — = te My (t) 


The following example looks at a case where X(s) is a rational function, 
but is not proper. 


Example 3.7 Find the inverse Laplace transform of 
Equation: 


X(s) s?+6s+1 
2 
s*+5s+6 


Here since gq = p = 2, we cannot perform a partial fraction expansion. First 
we must perform a long division, this leads to: 
Equation: 


- s—) 
s?+584+6 
s—) 
(s + 2)(s + 3) 


X(s): =1 


=1+ 


where s — 5 is the remainder resulting from the long division. The quotient 
of 1 is called a direct term. In general, the direct term corresponds to a 
polynomial in s. The partial fraction expansion is performed on the quotient 
term, which is always proper: 

Equation: 


s—5 Ay Ag 


(s+2)(s+3) st2. 843 


Using [link] gives 


Equation: 
Ae = s—9d 
s+ 3 s=—2 
==] 
Equation: 
s—9d 
a 
s+2 s=—3 
So we have 
Equation: 
C 8 
X(s)=1- 
(s) Bap 2 - s+3 
and 
Equation: 


Complex Conjugate Poles: 
Some poles occur in complex conjugate pairs as in the following example: 


Example 3.8 Find the output of a filter whose impulse response is 

h(t) = e~**u (t) and whose input is given by x(t) =cos (2t)u(t). Since 
the output is given by y(t) = x(t)*hA/(t), its Laplace transform is 

Y(s) = X(s)H(s). Therefore using the table of Laplace transform pairs 
we have 


Equation: 


S 
and 
Equation: 
1 
A = 
(s) s+5 
which leads to 
Equation: 
S 
2 
(s) (s? + 4)(s+5) 
7 S 
(s + 72)(s — j2)(s + 5) 
Aj A» Ag 


—f29) gaan Bee 


The poles are at s = 72, —j2 and -5, all of which are distinct, so equation 
[link] applies: 
Equation: 


A, =¥(s)(s +32) 


s=—]72 
_ s 
(s _ j2)(s = 5) s=—72 
__ 7 
—j4(5 — j2) 


5 + 72 
58 


The second coefficient is 
Equation: 


A, =Y(s)(s — j2)|,—j9 
5 — j2 
«+58 


The calculations for Az where omitted but it is easy to see that A» will be 
the complex conjugate of A, since all of the terms in A» are the complex 
conjugates of those in A;. Therefore, when there are a pair of complex 
conjugate poles, we need only calculate one of the two coefficients and the 
other will be its complex conjugate. The last coefficient corresponding to 
the pole at s = —5 is found using 

Equation: 


= 8 
(s? Flt 4) s=—5 
_ 8 
29 
This gives 
Equation: 
5+j2 5—j2 5 
( ) = 58 58 29 


We can now easily find the inverse Laplace transform of each individual 
term in the right-hand side of [link]: 
Equation: 


O92 «4 5 — 72 
_ 5 +22 a gy, 592, 


y (t) aa (t) = 


At this point, we are technically done, however the first two terms in y(t) 
are complex and also happen to be complex conjugates of each other. So we 
can simplify further by noting that 

Equation: 


SIZ 25 9-72 ; S92 4 
ey (t) + eu (t) =2Re (FP ety )) 
= 2Re(0. 0928e 1-38 67" (t)) 
= 0. 1857 cos (2¢ — 0. 3805) u(t) 


The simplified answer is given by 
Equation: 


y(t) = 0.1857 cos (2t — 0.3805)u (t) — 0.1724e~°*u (t) 


We note that the answer contains a transient term, —0. 1724e~‘u (t), and 
a steady-state term 0. 1857 cos (2¢ — 0. 3805). The steady-state term 
corresponds to the sinusoidal steady-state response of the filter (see Chapter 
3). It can be readily seen that the frequency response of the filter is 
Equation: 


1 
5+ 52 


H (j2) = 


and therefore |H(j2)| = 0. 1857 and ZH(j2) = —0. 3805. 


While the above example provides some insight into the sinusoidal steady- 
State response, the number of complex arithmetic calculations can be 
tedious. We repeat the example using an alternative expansion involving 
complex conjugate poles: 

Equation: 


1 Ajs+ Apo 


stbst+e s%+bst+e 


where it has been assumed that b? — 4c < 0 (otherwise, we have distinct or 
repeated real poles). As mentioned above, the expansion in [link] can be 
combined with expansions for distinct or repeated poles. 


Example 3.9 
Equation: 
8 
~ (s2 +4) (s +5) 
A;js+ Apo Az 
(s? + 4) s+5 


Using the cover up method gives 
Equation: 
9) 
an ol |e 
$2 + 4 s=—) 29 


Clearing fractions in [link] gives: 
Equation: 


s = (Ajs + Ap) (s +5) -— = (5? +4) 


Setting s = 0 in [link] gives Ag = , Substituting this value back into 
[link] and setting s = 1 leads to Ay = = The resulting Laplace transform 
is: 

Equation: 


5 4 
2ST a9 
(s? + 4) s+5 


Y(s) = 


Using the table of Laplace transforms then leads to 
Equation: 


Comparing this answer with [link], we see that the sum of a cosine and a 
sine having the same frequency is equal to a cosine at the same frequency 
having a certain phase shift and amplitude. In fact, it can be shown that 
Equation: 


a cos (Q,t) + b sin (Qt) = r cos (Qt — vy) 


with r = Va? + b? and (y =arctan 2. The following example also 


involves complex conjugate poles and illustrates some additional tricks to 
solving the partial fraction expansion. 


Example 3.10 Find the output of a filter whose input has Laplace transform 
X (s) = < and whose system function is given by 


Equation: 


Multiplying X(s) and H(s) gives 
Equation: 
1 
s(s? + 2s + 3) 
Ay Aos+ Ag 


Y(s) = 


Clearing fractions gives: 
Equation: 


1 =A, (s? + 2s + 3) +s (Ags + Az) 
= (A, + A»)s? + (2A, + A3)s+3A1 


Setting s = 0 leads to a quick solution for 4;, however two subsequent 
substitutions are needed to find Az and A3. A slightly faster way of solving 
for the coefficients in [link] is to rearrange the right hand side in terms of 
different powers of s (see second line). Then equate the coefficients of like 
powers of s on both sides of the equation to solve for the coefficients. For 
example equating the constant terms leads to 1 = 3.4, which gives 

Aya re The coefficients of s on either side of the equation are related by 


0 = 2A, + A3 which leads to A3 = ae Similarly, equating the 


coefficients of s? gives 0 = A; + Ay» which leads to Ay = — +. So we 
have: 
Equation: 
1 i 2 
Y(s) = 2 -—3- _3 


8 s* +2843 


The second term in Y(s) does not appear in most Laplace transform tables, 
however, we can complete the square of s* + 2s + 3 by taking one-half the 
coefficient of s, squaring it, then adding and subtracting it to give: 
Equation: 


e peel I= G41)" 42 


After a bit more massaging we get 
Equation: 


whose inverse Laplace transform is readily found from the table of Laplace 
transforms as 
Equation: 


y(t) = —u(t) - =e" cos (v2t)u (t) — e * sin (v2t)u (t) 


3/2 


Transfer Functions and Frequency Response 
Looks at the relationship between the transfer function and the frequency 
response. 


We saw in [link] that the transfer function of a linear time-invariant system 
is given by 
Equation: 


If we assume that H(s) is a rational function of s then we can write 
Equation: 


(s — 81) (s — Ba) --- (8 — By) 


Gea Genta) 


where (1, (2, ..., Bq are the zeros, and a1, Q2,..., @p are the poles of H(s). 
The poles and zeros are points in the s-plane where the transfer function is 
either non-existent (for a pole) or zero (for a zero). These points can be 
plotted in the s-plane with “x" representing the location of a pole and a “o" 
representing the location of a zero. Since 


Equation: 
|H(j2)| = H(s)| 30 
we have 
Equation: 
|g 2 — By||GQ — Bo|---|GQ— By| 
|H(g2Q)| = a 
|G2 — a4||GQ — ag|--- (GQ —ap| 
and 


Equation: 


Pp 


q 
LH(j2) = S~ Z(GQ_ — Bu) — Y| Z(GQR — 211) 


Each of the quantities 7{2 — 6; have the same magnitude and phase as a 
vector from the zero 6; to the j{2 axis in the complex s-plane. Likewise, 
the quantities 7{2 — a, have the same magnitude and phase as vectors from 
the pole a; to the 72 axis. 


Example 3.1 Consider a first-order lowpass filter with transfer function 
Equation: 


1 
A = 
(s) es 
Then 
Equation: 
H(j)| = — 
Ie Ti+ 2 


The quantity 7§2 + 2 has the same magnitude and phase as a vector from 
the pole at s = —2 to the 72 axis in the complex plane. The magnitude of 
the frequency response is the inverse of the magnitude of this vector. The 
length of 72 + 2 increases as 2 increases thereby making |H(j2)| 
decrease as one would expect of a lowpass filter. 


Example 3.2 Consider a first-order highpass filter with transfer function 
Equation: 


Then 
Equation: 


oy li 
HG) = 


When 2 = 0, |H(j2)| = 0, but as 2 approaches infinity, the two vectors 
9492 and j§2 + 2 have equal lengths, so the magnitude of the frequency 
response approaches unity. 


Example 3.3 Let's now look at the transfer function corresponding to a 
second-order filter 


Equation: 
1 
H(s) = ee ee 
(s+1+4 75)(s + 1-75) 
Or 
Equation: 
i2+1 
|H(J)| = sar 


\g72 + 1+ 75||G7QQ + 1 — 75] 


The pole-zero plot for this transfer function is shown in [link]. The 
corresponding magnitude and phase of the frequency response are shown in 
[link]. 


Pole-zero diagram for Example 
[link]. 


107 10°" 10° 10' 10° 
Frequency (rad/s) 
50 

or 
g oO 
o 
3 
8 
# -50 
a 

-100 

107 10° 10° 10' 10° 


Frequency (rad/s) 


Magnitude and phase of frequency response for Example 
[link]. 


The two poles are at s = —3 + 75 and the zero is at s = —1. When the 
frequency gets close to either one of the poles, the frequency response 
magnitude increases since the lengths of one of the vectors 72 + 75 is 
small. 


In the previous example we saw that as the poles get close to the 72 axis in 
the s-plane, the frequency response magnitude increases at frequencies that 
are close to the poles. In fact, when a pole is located directly on the 72 axis, 
the filter's frequency response magnitude becomes infinite, and will begin 
to oscillate. If a zero is located on the 72 axis, the filters frequency 
response magnitude will be zero. 


Example 3.4 Suppose a filter has transfer function 
Equation: 


s 


HAS) 55 


The two poles are at s = +75 and there is a zero at s = 0. From [link], the 
impulse response of the filter is h(t) =cos (5t)u(t). This impulse response 
doesn't die out like most of the impulse responses we've seen. Instead, it 
oscillates at a fixed frequency. The filter is called an oscillator. Oscillators 
are useful for generating high-frequency sinusoids used in wireless 
communications. 


Stability 
Looks at bounded input bounded output stability. 


Filter Stability 
We define a filter as being stable if a bounded input produces a bounded output. This is 


called BIBO stability. Consider the convolution integral 
Equation: 


We wish to find a condition that is necessary for the filter output y(t) to be bounded 
whenever the input z(t) is bounded. Let's take the absolute value of both sides of [link] 
Equation: 


ue)| =| fhe (ener 
< [ in@ete— nar 


< mas | |h(7)|dr < oO 


Co 


where the first inequality is due to the triangle inequality and the second inequality 
results from replacing |x2(¢ — 7)| by its upper bound 2 a2. Therefore, the condition for 
BIBO stability is that the impulse response the absolutely integrable. 

Equation: 


J ineoyat 2s 


Co 


Example 3.1 The impulse response of a filter is h(t) =cos (5t)u(t). If the input to this 
filter is x(t) =cos (5t)u(t), then the output is given by 
Equation: 


dey = i "eee (GR costo pues 


t 
/ cos (57) cos (5 (t — r))dr 
0 
Lr 1 ; 
= al cos (107 — 5t)dr + — cos (st) [ dr 
2 0 2 0 


Lf 1 
= sf cos (107 — 5t)dr + at cos (5t) 
0 


where in the third line, we have used the trigonometric identity 
Equation: 


cos (01) cos (82) = 5 (cos (8; — 02)+ cos (0; + 62)) 


As t approaches infinity, then it is clear that y(t) also becomes unbounded. Therefore 
the filter is unstable. Moreover, it is clear that h(t) is not absolutely integrable. 


If a transfer function is rational, it can be expressed as a sum of any direct term that 
may be present plus a proper rational function (q < p). 

Equation: 

bys? + by_18t | +--+ +b18 + by 


ApS? + Ap-18?-! +--+ +18 + ao 


1 
H (s) = cms” +Cm-18" +++: te1s +¢94 


The direct terms can be shown to produce unbounded outputs when the input is a step 
function (which of course, is bounded). The proper rational function produces output 
terms that depend on whether poles are distinct or repeated and whether these are real 
or complex: 


e Distinct real poles, s = ox lead to an impulse response with terms: 
Equation: 


Ke u(t) 


¢ Distinct complex conjugate poles s = o% + 7§2, produce impulse response terms: 
Equation: 


Ke™ cos (Qyt + 0)u(t) 


e Repeated real poles s = o, produce impulse response terms having the form: 
Equation: 


Kt"e™*u (t) 


The quantities K and @ are constants while n is a positive integer. In all of these cases 
the filter impulse response dies out with time if the poles have negative real parts, i.e. it 
is absolutely integrable and therefore leads to a stable filter. If the poles have zero or 
positive real parts, then the impulse response terms either oscillate or grow with time, 
and are not absolutely integrable. When this happens then one can always find a 
bounded input that produces an unbounded output. Therefore, in order for a filter with 
a rational transfer function to have BIBO stability, the transfer function should be 
proper and the poles of the transfer function should have negative real parts. 


Second-Order Filters 
Summarizes second order lowpass, bandpass, and highpass fitlers. 


Second-Order Lowpass Filters 


The second-order lowpass filter has system function 
Equation: 


Q2 


2G) == 
(8) = yon, + 


where £2, is the undamped natural frequency, and ¢ is the damping ratio. 
The undamped natural frequency and damping ratio are properties of the 
physical devices used to implement the second-order filter (capacitors, 
inductors, and resistors, in the case of an electrical circuit). It so happens 
that the damping ratio satisfies ¢ > 0. Using the formula for the roots of a 
quadratic polynomial, the two poles of H(s) are easily found to be 
Equation: 


81 =—C2%m+ MVC? —-1 
8. =—C€Qn— MVC? -1 


There are three possible sets of poles that are categorized as follows: 


1. Overdamped: the poles are real and distinct. This occurs if ¢ > 1. In 
this case the impulse response is given by: 


Equation: 
(2, i254 
h(t) = ———e*""u (t) - ———e** u(t) 
2,/C2 —1 2) Coal 


2. Critically damped: corresponds to ¢ = 1. The two poles are repeated 
with, 
Equation: 


S| — 42.—> =C12,, = =D 


and the impulse response is given by 
Equation: 


h(t) = 27 te" (t) 


3. Underdamped: corresponds to 0 < ¢ < 1, giving a pair of complex 
conjugate poles. In this case, the impulse response is given by 
Equation: 


nin) = Pn cant sin (V1— Cnt) u (t) 


ae 


Note that in the underdamped case, the magnitude of the poles is 

|$1 |=|$2|= 2,, and when ¢ = 0, the two poles are on the imaginary 
axis which corresponds to an impulse response that is a pure sinusoid 
and the system is unstable. 


A root locus diagram shows the paths that the poles of H(s) would take as 
the damping ratio is decreased from some number greater than 1 down to 0, 
and is shown in [link]. 
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Root locus of second-order lowpass filter 
having 2, = 10. 


The frequency response of the second-order lowpass filter can be found 
using the substitution H (j2) = H (s)|,_;9, giving 


Equation: 

Q?2 

HGQ) = ——_"*» 
U2) = oy 924 pace, N 


The frequency response magnitude is shown in [link] for 2, = 10 and 
several values of ¢. Note that for the underdamped case, there is a 


resonance or peak that would be expected to maximize at the undamped 
natural frequency (2,,, when ¢ = 0. 


10° 10° 10' 10° 


C2 


Frequency response of second-order lowpass filter with 22, = 10 
and several values of ¢. 


We also observe that for the critically damped case, since 
Equation: 


Q?2 
(8+ Qn)” 


setting s = j gives |H(2,)| = +. Moreover, it is clear that for 2 < Qn, 
|H(2,,)| > + and for 2 > Q,, |H(Q2,)| < +. These ideas will be used 
later when we discuss Bode plots. 


Second-Order Filter Implementation 


Series RLC circuit. 


Second order filters can be implemented using either passive or active 
circuit elements. We will consider the series RLC circuit shown in [link]. A 
lowpass filter results by taking the filter output to be the voltage across the 
capacitor. Using voltage division, the system function is easily found to be 
Equation: 


Vip (s) 


Hip (s) = 


Comparing [link] with [link] we find that the undamped natural frequency 


is (2,.= a while 2¢2, = a giving the attenuation coefficient: 


Equation: 


A second-order highpass filter is implemented by taking the output of the 
filter to be the inductor voltage in the series RLC circuit. The resulting 
system function is 

Equation: 


Hip (8) = 7 


So in terms of the undamped natural frequency and damping ratio, the 
second-order highpass filter is given by 
Equation: 


s2 


Hi (8) = 3c + 


The graph of the frequency response magnitude of this filter is shown in 
[link] for several damping ratios and (2,, = 10. 


|H(j2)| (dB) 


Frequency response of second-order highpass filter 
with §2,, = 10 and several values of ¢. 


Not surprisingly, a second-order bandpass filter results by taking the output 
of the filter be the resistor voltage in the series RLC circuit. The system 


function is then 
Equation: 


Hyp(s) = 2 


This can be expressed in terms of ¢ and §2,, as 
Equation: 


262y8 


aera oy 


Frequency Response of Second-Order Bandpass Filter 


Setting s = 7/2 in [link] gives the frequency response of the second-order 
bandpass filter 
Equation: 


J262, 2 
Hee GO) = oa 1+ 7262, 


The magnitude of this frequency response is shown in [link] where 
7. = 1). 


| H(j2)|(aB) 
8 


-50 


Q 


Frequency response of second-order bandpass 
filter with (2, = 10 and several values of ¢. 


Evidently, the frequency response peaks at {2 = §2,. To prove this, we can 
divide the numerator and denominator of [link] by 72¢2,,22, giving 
Equation: 


Hyp (j22) — 


where 
Equation: 


2? — 02 


A= Sea 


It is easy to check that | A(2,,)| = 0 and |A(2)| > 0, 2 ~ 2,,. Therefore 
| Hop (72)| does in fact peak at Q = Qy. 


The bandwidth of bandpass filters is typically measured as the difference 
between the two frequencies, {22 and §2;, at which the gain has dropped by 
3 decibels from its peak of 0 dB. A diagram of this is shown in [link]. 
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Definition of bandpass filter bandwidth, 
BW = 22g — 921. 


Formulas for (2; and {22 can be found by solving 
Equation: 


since if A(2) = +1, we get |H(j2)| = 


A(Q) = +1 


1 : ‘ 
Ts: Solving the quadratic 


equation A({2) = 1 gives two solutions: 
Equation: 


Equation: 


Qo = CA — Qn V1 + CP 


while solving A(2) = —1 gives two additional solutions 
Equation: 


3 = —62y, ate QnvV/1 a Ce 


Equation: 


24 ms =¢12,, - QnV1 = ce 


Since 22,, af: 1+ ¢2 > CQn, we pick the two positive solutions as the 
“corner” frequencies: 
Equation: 


2, = —CO, + AVI + CP 


Equation: 


Qy = CQ, + MmVI+C 


The two negative frequencies are just the negatives of the two positive 
frequencies (recall that |H(j2)| has even symmetry). So the bandwidth is 
BW = 22 — 2) = 2CM2,. It can also be readily verified that 

Qn = V21Qo, that is, the center frequency of the bandpass filter is the 
geometric mean of the two corner frequencies. The qualify factor, Qo is 
defined as 

Equation: 


and is a measure of the narrowness of the filter bandwidth, with respect to 
its center frequency. It can be seen that Q, = ae 


Bode Plots 
Describes how to manually sketch Bode amplitude and phase plots. 


Bode Plots 


Plotting the magnitude and phase of the frequency response is most easily 
accomplished with a computer, provided you have the right software (for 
example, the Matlab function “freqs”). However if there is no computer handy, 
a classical method for quickly sketching the magnitude or phase response of a 
filter is using a Bode plot. Consider a general system function given by 
Equation: 


a (s pa 


H(s = ky —————_. 
a (s — a)"* 


where we assume the ‘yz and 6, are integers. The corresponding frequency 
response is therefore 


Equation: 
| pat (J82 — Be)” 
H (jQ) = eee 
il (72 _ az) 


[link] can be written as 
Equation: 


where 
Equation: 


4 By (-1)™ 


5 5 
= aad)" 


k=k, 


Since most frequency response plots are expressed in units of decibels, we 
have 
Equation: 


20 logiy |H(G2)| = 20 log, K 


Jk 
a?) 
kel 1-j2) 
= 20 logy |A| 
q Q Pp Q 
+S 7420 logig 1- iB. — S | 6420 log;y, 1- Iz 
k=1 k=1 


where we have used basic properties of logarithms. The individual terms in the 
sums can be plotted, to within a reasonable approximation, with relative ease. 
So the Bode magnitude plot involves summing the graphs of each individual 
term in [link]. Lets consider first a single positive term having the form 
Equation: 


2 
M = 720 log;, 1- IB 


It is clear that when 2 < 6 then M ~ 0. On the other hand if 2 > £, 

M = 720 logig 3 = 720 log, |Q| — 720 logy, |G]. If we plot this as a 
function of log,, |{2|, this represents a straight line having a slope of y20 that 
crosses the logy, |2|-axis at logy, |§2| =log,, |G]. These approximations are 
illustrated in [link]. 
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Straight-line approximations to 


Instead of plotting M as a function of logy, ||, it is more common to plot it 
as a function of (2 with the frequency axis on a logarithmic scale. In this case 
the non-zero slope is 207 decibels per decade, where one decade represents an 
increase in {2 by a factor of 10. The modified graph is shown in [link]. 
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Straight-line approximations to 
720 log;, 1—-— ig using logarithmic 
frequency axis. 


log scale 


We also note that when {2 = 8, the straight-line approximations are not valid, 
however the true value is easily found to be 


720 log,, |1 — j|/ = 720 log,, V2 = 34 dB. Negative terms in [link] are 
approximated in a similar manner, but the non-zero slope is now —207 
dB/decade. The resulting approximation is shown in [link]. 
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Straight-line approximations to 
—§20 logy, |1— j4|. 


[link] did not take into account cases where the poles or zeros occur at s = 0. 
For example, if H (s) = s7, then M = 20 log), |AH(jQ)| = 720 logy) |Q, 
which is a line having a slope of 20y dB/decade passing through the 2-axis at 
§2 = 1 (see [link]). When there is a single or repeated pole at the origin, the 
graph appears just as in [link] but with a negative slope. 
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Straight-line approximations to 720 log,) ||. 


Next we'll look at approximations to the phase response. Here we'll begin with 
H(jQ2) as shown in [link]. Taking the phase of both sides gives 
Equation: 


Mea (1 - i) 
ie) 
= 2K + c(1-72) - Yo hZ (1-i=) 


k=1 


ZH(jQ) =Z K 


where we have used the fact 7{Z1Z2} = 2Z, + ZZq and Z{Z7} = yZZ. 
Lets now consider how we approximate each of the terms in [link]. Consider 


the single term 
2 2 
Z\1 -35 = — arctan @ 
( B B 


Equation: 
The graph of this function is shown in [link]. 


Graph of — arctan (4) , shown on a linear frequency 


scale. 


Since the magnitude response plots are on a logarithmic frequency axis, it 
would be desirable to do the same for the phase response plots. If we restrict (2 
to positive values and plot — arctan 


i ) on a logarithmic frequency scale, 
we get the graph shown in [link]. 


0 = 


Q 


Graph of — arctan (4) , shown on a logarithmic 
frequency scale, for (22 > 0. 


This graph can be approximated with straight lines as shown in [link]. The 
approximation is a straight line having a slope of —7/4 rad/decade passing 
through a phase of —7/4 at 2 = £. The line then levels off at (2 = 0.16 and 
§2 = 108. The resulting approximation is shown in [link]. 


0.13 B 103 


Approximation (in red) to graph of — arctan (4) ona 


logarithmic frequency scale. 


Constant terms in [link] have a phase of 0 or +7, depending on their sign, 
while poles or zeros of H(s) at the origin produce phase terms of 7/2 for 
zeros or order ‘y or —d7/2 for poles of order 6. Next we'll illustrate these 
techniques with a few examples. 


Example 3.1 Sketch the Bode magnitude and phase response plot for the 
following filter: 
Equation: 


7 
(s + 10) (s + 10°) 


We begin by finding the corresponding frequency response by setting s = 7/2: 
Equation: 


i 1 
H (jQ) = (+8) (0458) 


Lets find the magnitude response first: 
Equation: 


Q2 Q2 


The resulting straight-line approximations to the two terms in [link] along with 
their sum are shown in [link]. 
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Bode magnitude response plot derivation for 
Example 1. 


The phase is given by 
Equation: 


i 2 


The straight-line approximations to the two phase terms, along with their sum 
are shown in [link]. 
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Bode phase response plot derivation for Example 1. 


Next we'll look at an example of a second-order highpass filter. 


Example 3.2 Find the Bode magnitude and phase response of the following 
filter 


Equation: 
2 
H(s)= Ee . 

(s + 10) 
Substituting s = 72 in [link] gives 
Equation: 

—10-°2? 
H (G2) = ~——_ 
(1+ 49) 


The magnitude response, expressed in decibels becomes 
Equation: 


2 
20 logig |H(j2)| = 20 log,, 10-2 +40 logy, |2| — 40 log,, 1+ a 


The graphs of the straight-line approximations for the three terms in the right- 
hand side of [link] along with their sum are shown in [Link]. 
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Bode magnitude response plot derivation for 
Example 3.2. 


The phase of the frequency response is found to be 
Equation: 


LH(jQ) = Z(-10-72?) — Z (1 + 


The resulting Bode phase response plot is found in [link]. 


2H (jQ) = 7-22 (1+ 2) 
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Bode phase response plot derivation for 
Example 3.2. 


The Laplace Transform: Excercises 


Exercises 


1. Find the Laplace Transform of the following signals, for each case 
indicate the Laplace transform property that was used: 


aot) ae a (2) 

b. x (t) = 4te~°"u (t) 

c.  (t) = 4e~9 24-1), (¢ — 10) 
d. x(t) = d(t — 5) 

e. z(t) = 10tu(t) 

f. x(t) =sin (10zt)u(t) 

g.2(t) =e * sin (107t)u (t) 
h. a(t) = rect(t — 0.5, 1) 


2. Suppose that two filters having impulse responses h, (t) and hg (t) are 
cascaded (i.e. connected in series). Find the transfer function of the 
equivalent filter assuming h, (t) = 10e~!'u (¢) and 
hy (t) = Seu (t). 

3. Find the inverse Laplace transforms of the following: 


s)= +> 


(s) = 
(s) = 35 
( = 1 


a. X 

beaX. 

c. X (a) (s+3)? 
d. X(s) = 10 

e. X(s)= 4 
ee Se) e— 


== 


4. Use partial fraction expansions to find the inverse Laplace transforms 
of the following: 


-_ s+2 
a. X (8) = Geyser) 
a” X(s) = s+l1 


(s+2)°(s-+3) 


10. 


 X(s) = Geary 


——_-g?—3s+1 
d. X (s) = (e643) 


. Consider a filter having impulse response h (t) = e *“u (t). Use 


Laplace transforms to find the output of the filter when the input is 
given by: 


. Indicate whether the following impulse responses correspond to stable 


or unstable filters: 


a. h(t) = u(t) 

bh(tiS er ae) 

c. h(t) = e~* cos (4t)u (t) 
d. h(t) =cos (10t)u(t) 


. Use Laplace transform tables to find the impulse response of the 


second-order lowpass filter in terms of ¢ and 2, for the overdamped, 
critically damped, and underdamped case. 


. Use a series RLC circuit to design a critically damped second-order 


lowpass filter with a corner frequency of 100 rad/sec. Usea R = 6.8 k 
{2 resistor in your design. 


. Using a 10 k2 resistor, design a critically damped bandpass filter, 


having a center frequency of 100 rad/sec and indicate the resulting 
bandwidth of the filter. What is the quality factor of the filter? 

Use bode plots to find the magnitude and phase response of the 
following filters 


a. Equation: 


b. Equation: 
10°s 


ao (s + 10)(s + 10°) 
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